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Chapter 1 
Introduction 



Carbon is a basic ingredient of life and all organic chemistry which is con- 
sequence of its abundance in nature and his chemical reactivity. With four 
valence electrons distributed to one 2s and three 2p orbitals, which can hy- 
bridize in many different ways, carbon is characterized by a large flexibility 
of chemical bonding. One particularly interesting case is sp2 hybridization 
which creates three strong cr-bonds in plane, while the remaining p orbital is 
weakly bound with neighboring atoms creating 7r-bond. In this thesis we will 
be studying graphene: a two-dimensional (2D) crystal of carbon atoms assem- 
bled in a honeycomb structure. While cx-bond is responsible for the most of 
the structural integrity of graphene, 7r-bond determines low-energy electric and 
optical properties. Very peculiar property of graphene is that its low-energy 
electrons behave as massless Dirac particles [TJ [2] (near the corners of the Bril- 
louin zone). Since graphene is essentially a 2D material, one can simply tune 
its Fermi level through an electrostatic gating which brings about large con- 
trol over electrical and optical properties, important for various technological 
applications. 



1.1 Experimental realization 

Scientists were puzzled for long time whether nature allows existence of a two- 
dimensional crystal. In 1930's Peierls [3] and Landau [I] showed that thermal 
fluctuations would destroy long range order and essentially melt 2D lattice 
at any finite temperature. Therefore it came as a surprise when Geim and 
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Novoselov announced El [7] in 2004 a discovery of a first 2D crystal made 
of carbon atoms - graphene. Scientist were further astonished by a shear sim- 
plicity of the experimental method which essentially used a scotch tape to 
exfoliate graphite (graphite can be viewed as a simple stack of weakly bound 
graphene planes). The 2010 Nobel prize in physics came as a credit for this 
great discovery but it is interesting that even today in 2012 experimentalists 
still use this " scotch tape technique" since it offers exceptionally pure graphene 
samples on a small scale, important for fundamental research. Of course it is 
impractical on a large scale production which is required by various industrial 
applications, and soon after the discovery of graphene several other methods 
were developed for graphene production, most notably chemical vapor deposi- 
tion (CVD) [H], segregation by heat treatment of carbon-containing substrates 
[9] and liquid phase exfoliation [10] . The most promising of these methods, for 
large scale graphene growth, is CVD which is also used [11] by the group of 
Dr. Marko Kralj from the Institute of Physics in Zagreb, Croatia. They heat 
ethylene (C2H4) gas, up to a temperature of 1000°C, above the metal surface 
which serves both as a catalyst for ethylene decomposition and substrate for 
graphene growth. 

It is interesting to note that various groups claim they have seen graphene in 
their experiments prior to 2004 but it wasn't until Geim and Novoselov ground- 
breaking experiments that the true potential and importance of graphene was 
recognized. 

While graphene's intriguing mechanical properties are still debated, this the- 
sis concerns primarily with electrical and optical properties which are a subject 
of intense research and numerous practical applications. 

1.2 Plasmonics 

Plasmonics studies collective electron surface charge oscillations (surface plas- 
mons at surfaces of bulk materials or plasmons in a pure 2D materials like 
graphene) accompanied by tight confinement of electromagnetic (EM) fields. 
In recent years, an enormous interest has been surrounding the field of plas- 
monics, because of the variety of tremendously exciting and novel phenomena 
it could enable. On one hand, plasmonics seems to be the only viable path 
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toward realization of nanophotonics: control of light at scales substantially 
smaller than the wavelength [121 Q21 US] . On the other hand, plasmonics is 
a crucial ingredient for implementation of most metamaterials, and thereby all 
the exciting phenomena that they support [16], [T71 [HI [19] , including negative 
refraction, superlensing, and cloaking. However, there is one large and so far 
insurmountable obstacle towards achieving this great vision: plasmonic mate- 
rials (most notably metals) have enormous losses in the frequency regimes of 
interest. This greatly motivates us to explore plasmons and their losses in a 
newly available material with unique properties: graphene [SI El E] - 

Plasmons are also very interesting phenomenon from the point of view of 
many-body physics. Since losses are in a large manner determined by phonons 
we will encounter interactions between various elementary excitations and in- 
teresting many-body effects like plasmon-phonon coupling. 

1.3 Near field thermo-photo-voltaics 

Radiative heat transfer between two bodies can be greatly enhanced in the 
near field, i.e. by bringing the surfaces close together to allow tunneling of 
evanescent photon modes [201 ED E2] • This happens because near field radi- 
ation transfer involves thermal excitation of various surface modes which can 
have much greater wave vectors (and density of states) than the freely propa- 
gating modes (limited by the light line). Since each wave vector corresponds 
to a heat channel, vacuum becomes better heat conductor in the near field. 
However, due to their localization and evanescent nature, it is only at sub- 
wavelength separations that these modes become relevant. While measuring 
near field transfer has been experimentally difficult [231 EH EH [26] , the promise 
of order-of-magnitude enhancements over the far field Stefan-Boltzman black 
body limit has made transfer in the near field the topic of much research. 

With the current world energy demand and large environmental impact of 
fossil fuels there is a worldwide shift toward renewable energy sources. In 
that respect, thermo-photo-voltaics (TPVs) are a promising class of heat to 
electricity conversion devices [271 EH] where Sun can heat up an emitter that 
selectively re-radiates frequencies matched to the band gap of the photo- voltaic 
cell thus minimizing the thermalization losses. TPVs are not limited by the 
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Sun source and can use any hot (terrestrial) object like a factory furnace or 
various hot car parts as a heat source. From the perspective of future energy 
crisis there is a large demand for more efficient energy management where 
TPVs can play important role by turning wasted heat into electricity. 

Near field TPVs j29l (30j EI] further offer greater power densities since the 
near field heat transfer can be orders of magnitude larger than the far field 
limit. Finally, due to evanescent nature of EM modes, one does not need to 
worry about losing energy through modes with frequencies below the photo- 
voltaic band gap, resulting in even larger device efficiencies. 

1.4 Objectives and results 

The objective of this research is to study electrodynamic properties of graphene 
and especially high-frequency collective oscillations of electrons (plasmons). 
We will analyze plasmon excitations in the context of plasmonics and other 
technological applications, but we will also look at the same problem from 
the point of view of many-body physics as an interaction between various 
elementary excitations (plasmons, phonons, etc.). Finally we study near field 
heat transfer with graphene (mediated by thermally excited plasmons) in the 
context of TPVs. 

We study plasmon excitations in graphene in the context of the Random 
Phase Approximation (RPA) [56] and number-conserving relaxation-time ap- 
proximation [34J and we show that plasmons in doped graphene can have both 
low losses and large localization for frequencies below optical phonon energy at 
0.2 eV. Large plasmon damping occurs in the regime of interband single parti- 
cle excitations which can be shifted towards larger energies for stronger doping 
values. We demonstrate that for sufficiently large doping there is a frequency 
interval from optical phonon frequency to boundary of interband regime, where 
the plasmon damping is dominated by emission of optical phonon and electron- 
hole pair. To describe impurity scattering we use DC relaxation time since we 
don't expect significant frequency dependance. The phonon contribution is es- 
timated from the electron self-energy induced by electron-phonon interaction. 

We also explore electron-phonon interaction in graphene as an interesting 
problem from the aspect of many-body physics. By measuring Raman shift of 
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optical phonon energy it was demonstrated that Born-Oppenheimer approxi- 
mation (BOA) is not a valid approximation in graphene [32J. The measured 
Raman shift is a consequence of the interaction with single particle excitations, 
however the breakdown of BOA means that electrons and phonons move on 
comparable energy scales which leads to a possibility of interaction between 
phonons and collective electron excitations (plasmons). We show that a pecu- 
liar type of hybridization of plasmon and optical phonon modes occurs around 
the point where the two modes cross in energy and momentum simultaneously 
since then the electron-phonon interaction will be drastically increased due to 
collective electron response. We demonstrate that the electron-phonon interac- 
tion leads to polarization mixing of the two modes so that longitudinal plasmon 
(LP) couples exclusively to the transverse optical phonon (TO) mode, while 
the tranverse electric mode, also referred to as the transverse plasmon (TP), 
couples exclusively to longitudinal optical phonon (LO) mode; thus there is no 
coupling between LPs and LO modes. Formally, we analyze plasmon-phonon 
coupling in the self-consistent linear-response formalism which describes inter- 
action of phonons with both single particle and collective electronic excitations. 
We emphasize that the phonon interaction with collective excitations is much 
larger than the phonon interaction with single particle excitations (measured 
by Raman) which means that plasmon-phonon interaction can serve as a mag- 
nifier for exploring electron-phonon interaction in graphene. Further on, our 
calculations give a slight correction to the standard result of Raman shift of 
the optical phonon energy since the longwave phonons can interact also with 
radiative EM modes so that we predict increasing Raman linewidths for higher 
dopings. Finally we note that LO phonon decouples from all (single particle 
and collective) electronic excitations when its dispersion crosses the light line. 

While longitudinal charge density oscillation can be referred to as longitudi- 
nal plasmon, which is also polarized like transverse-magnetic (TM) EM mode, 
we also analyze properties of the unusual transverse plasmon in 2D systems 
|48j . which is polarized like transverse-electric (TE) mode, and accompanied 
by transverse current density oscillation. These kind of modes are possible 
only if the imaginary part of 2D conductivity is negative which in principle 
requires interband transitions. From that perspective bilayer graphene is an 
interesting candidate for exploring these modes, because it has a rich band 
structure and particularly two perfectly nested bands with a gap of 0.4 eV 



5 



which results in large joint density of states considering the vertical interband 
transitions. We show that plasmon properties (localization) of TE modes are 
much more pronounced in bilayer than in single layer graphene. 

We also show that thermally excited plasmons strongly mediate and enhance 
the near field radiation transfer between two closely separated graphene sheets. 
Near field heat transfer is analyzed within the framework of fluctuational elec- 
trodynamics and we predict several orders of magnitude larger values of heat 
transfer between two graphene sheets in the near field than the case of heat 
transfer between two black bodies, of the same temperatures, in the far field. 
Finally we demonstrate that graphene can be used as a thermal emitter in the 
near field thermophotovoltaics leading to large efficiencies and power densities. 

The thesis is organized into chapters as follows. In Chapter [2] we present 
theoretical methods and tools that will be used throughout the text. We first 
calculate electron dispersion and electron-phonon interaction Hamiltonian in 
graphene within the tight binding approximation. Next we give the density- 
density and current-current response functions in the linear approximation 
and use fluctuation-dissipation theorem to calculate current- current correlation 
function due to thermal fluctuations in the system. Finally we use this to 
calculate the radiative heat transfer between two graphene sheets. In Chapter 
|3]we calculate plasmon dispersion and damping due to electron-impurity and 
electron-phonon scattering. In Chapter H] we calculate dispersion od TE modes 
in single and bilayer graphene. In Chapter we calculate plasmon-phonon 
interaction within the self-consistent linear response formalism. In Chapter 
|6] we calculate near field heat transfer between two graphene sheets and we 
analyze near field TPV device with graphene as a thermal emitter. Finally, in 
Chapter [7| we summarize. 
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Chapter 2 
Methods 



In this chapter, for the sake of the clarity of the presentation, we derive basic 
physical quantities used to describe graphene such as the low energy Dirac 
Hamiltonian and the electron-phonon interaction. We will also define standard 
response functions, like the conductivity, density-density, and current-current 
response functions, that will be used in later chapters. This chapter is intended 
to provide an introduction and overview of these concepts so the reader already 
familiar with them can skip the corresponding sections. Finally we will derive 
an expression for the radiative heat transfer between two graphene sheets at 
different temperatures by employing the fluctuation-dissipation theorem. 

2.1 Tight binding approximation in graphene 

In this section we use the tight-binding approximation to derive the electron 
band structure of graphene, Dirac equation valid at low energies and electron- 
phonon interaction. 

2.1.1 Electron band structure 

Graphene crystal structure is determined by a Bravais lattice with two atoms 
in a basis (see figure [XT]) . We can choose unit cell vectors as sli = a(l,0) and 
a 2 = a (—1/2, \/3/2), while the vectors connecting first neighbors are given by 
Ti = a(0,l/V3), T2 = a(— 1/2, —1/2-^/3), and T3 = a(l/2, — l/2y/3). Here 
a = 0.25 nm is a lattice constant while the nearest neighbor carbon-carbon 
distance is |t/| =b = a/\/3 = 0.14 nm. 
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Unit cell vectors of reciprocal lattice are given by Ci = (2n/a)(l, 1/y/S) 
and c 2 = (27r/a)(0,2/\/3), while we are primarily interested in the vertex 
points of the Brillouin zone i.e. vectors K = (2?r/a)(l/3, 1/V3) and K' = 
(27r/a)(2/3, 0). The remaining four vertex points are equivalent to the points 
K i K' since they are connected to them by a simple translation with the 
reciprocal vector n\C\ + ?^ 2 c 2 , where n\ and n 2 are integers. 



(a) (b) 




Figure 2.1: a) Graphene crystal structure. Unit cell vectors are ai and &2 
while A and B are atoms of the basis, b) Brillouin zone. We mark high sym- 
metry points K and K' where the low-energy electron excitations are described 
by massless Dirac equation. 

As we already pointed out in the introduction, sp2 hybridization is responsi- 
ble for the mechanical stability of graphene by creating three strong a bonds in 
xy plane while the remaining p z orbital weakly interacts with the neighborings 
p z orbitals creating the ir bond. Since we are particularly interested in tt bond, 
the entire problem is very well described with the tight binding approximation 

PP- 

Let us define now an operator that creates a free p z orbital at the lattice 
point R, i.e. \p z (TVj) = CjjO). Let us further denote by —70 the hopping inte- 
gral between nearest neighbor p z orbitals (next-nearest neighbor interaction is 
negligible and 70 ~ 2.8 eV [I]). Since we are only interested in the behavior of 
the electron energies near the p z orbital energy, our system is well described 
by a tight binding Hamiltonian 

#=-70 CRa-T^a-To 4 b +T i R b > t 2 ' 1 ) 

Ra,T, Rb,T, 

where the sum over lattice points is divided into two parts that contain different 
basis atoms i.e. Ra = ^1^1 + n 2 a 2 + T\ and Rb = n\B.\ + n 2 a 2 {ri\ and n 2 are 
integers). In equation (12. ip . we have assumed that zero energy corresponds to 
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p z orbital energy (i.e. E(p z ) = 0) and we have neglected overlapping of the 
two neighboring orbitals. We have also omitted the notion of electron spin 
since it only plays the role of additional degree of freedom. Eigenstates of the 
Hamiltonian (12. ip must take the form of the linear combination of p z orbitals 
that satisfy the Bloch condition 

4 = ^ E ^ RA ck A ^(k) + e^dj B (k)Z, (2.2) 

Ra Rb 

where we have explicitly separated the phase Z (Z*Z = 1) which will be 
defined later so that analytical expressions would look as simple as possible. 
Let us define now Fourier transform of operators c^ A and c^ B as 

4 = 4^^4^ ( 2 - 3 ) 



^ = ^£ elkRB 4 B - (2-4) 

V Rb 

Then, the Bloch eigenstate Q is c[ = f A (k)A j k + f B (k)Bl, and we can also 
write the inverse Fourier transforms (since Z*Z = 1) as 



Cra= vn 



J=5>- k - R -4,and (2.5) 



Z* 

c 



Rb 

■ 



Let us look now at the first sum from (12.11) and notice that every vector Ra — t~i 
is in fact one of the Rb vectors, so we have 

£ <L-T,c** = £ (^£e-— »^) (^£~ 

= Z*e ik ' T 'BlA k J2^ el{k ^' )RA - ( 2 - 7 ) 

k,k',T i R A 

However, since 

E^ ei(k " k ' )RA = ^, (2.8) 

R A 
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we obtain in the first sum 

E 4 A -r ( % A = E Z*^ T 'BlA k . (2.9) 

Ra,T, k,T( 

In a similar manner we get the second sum 

= E Ze-^Aj^k. (2.10) 
Finally, the Hamiltonian (12. ip becomes 

H = -to E ( E z * eikTi • ^ + E Ze_ik ' Ti • ) • (2.H) 

k V i i j 

Relation (12. lip contains a specially important function 

T(k) = -7o E Ze ~ Jk ' Ti > ( 2 - 12 ) 
so finally we can write equation (12. lip in a matrix form 

--?(4*)U7)(i)- 

Now, since the Bloch state = /^(k)^ + /^(k)5^ has to diagonalize this 
Hamiltonian, we can also write 

h = y,e^cW = X>w (4 * ) ( JJJ ) ( am ft« ){ A b\ 

(2.14) 

By comparing equations (12.131) and (12.141) we need to have: 

( ° ^ )( = B(k) ( V (2.15) 

^ T-(k) / ^ / B (k) ^ \ / B (k) / 

So we have reduced entire problem to the matrix diagonalization, while the 
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Figure 2.2: Graphene electron band structure with Dirac cones around K 
point (magnified). Intrinsic graphene has Fermi level Ep = = 0. 



eigenvalues (i.e. energies) are given by: 

= 0. (2.16) 



-E(k) T(k) 
T*(k) -E(k) 



Solution of the determinant equation (I2.16P determines the electron band struc- 
ture in graphene as [I]: 



E±(k) = ±v/T(k)T*(k) = ±70^ 1 + 4 cos ^ cos + 4 cos* ^. 

(2.17) 

Figure 12721 shows the function E(k) and we can notice the peculiar behavior 
of the bands at the Brillouin zone vertex points K and K'. Further on, since 
each graphene unit cell contains two atoms in basis and each atom donates one 
free electron into the band, Fermi energy is defined such that there are enough 
electrons to fill precisely one Brillouin zone in the reciprocal space. Relation 
(12.171) tells us that electron bands are divided into positive and negative states 
that touch precisely at the the vertex point of the Brillouin zone (see also 



figure 1272]) . such that we have Ep = -Ek = 0. Because of the fact that electron 
states around the Fermi energy determines the low-energy properties, we will 
focus precisely on the area around the K and the K 7 points. Finally we note 
that, since valence (negative) and conduction (positive) band touch at only 
6 points (K, K' and the remaining four equivalent vertex points), that are 
located precisely at the Fermi level, the intrinsic graphene is an unusual zero 
gap semiconductor. 



11 



2.1.2 Dirac electron dispersion in graphene 



We can write equation ( I2.15P as an eigenvalue equation: H-^ip-^ = E k ip k , where 
the Hamiltonian and the wave function (eigenfunction) are given by 

H k =( ° T(k) Vand (2.18) 
\ T*(k) J 

Let us focus now on the area around the K point and change the origin of our 
wave vector as k — y k + K, so that we have |k| << |K|. Now we can make a 
Taylor expansion of the function T(k) as follows: 

T(k) = -7o Ze~ t{k+K) - Tl « -7o ^ Ze- l ^ Tl {\ - ik ■ r t ). (2.20) 

i i 

Next we calculate the following sums: 

e~ iK - T ' = 0, and (2.21) 



]T ne-^ Tl = e -^ 3 ^-a(tx + y). (2.22) 
i 

Now we will choose the phase Z = e~ i7r / 3 so that we have 

T L Ze-^ = e-^a(ix + y) = ^fa(-ix - y). (2.23) 



Finally we get expressions for the function T(k) and the effective Hamiltonian 
in the vicinity of point K: 

/3 

T(k) = 7o zk • riZe~ iK - Tl = ^-a lo (k x - ik y ), and (2.24) 



/ k x — ik v 

* = V»U + tt. o (2 ' 25) 
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It is now convenient to introduce new variable: hvp = ^070, where vp ~ 10 e 
m/s since 70 ~ 2.8 eV pp. Hamiltonian f !2.25|) now becomes 



hvi 




(2.26) 



that is, 



where cr = cr^x+cryy, while a x 



= hvF<r ■ k, 
1 



and a,, 







(2.27) 



are the Pauli 



1 J \ % 

spin matrices. Here we note the remarkable property of graphene around K 

point where electrons behave precisely like massless Dirac particles of spin 1/2 

[53] ! We can also find energies (eigenvalues) and wave functions (eigenvectors) 

from the equation H^ip^ = E^tp^: 



E n ,k = n ■ Hvp\k\ = n ■ YwF\jk% + ky,and 



(2.28) 



VVk(r) = (r|nk) 




ik-r 



(2.29) 



Here L 2 is the area of graphene, n = 1 (n = —1) denotes the conduction (va- 
lence) band, respectively, and the angle 8(k) = ta,n~ l (k y /h x ). Further on, we 
note that behavior around K' point is easily found if we move the wave vector 
origin so that k — y k + K'. In that case it is more convenient to choose the 
phase Z — 1 and the Hamiltonian f 1 2 . X 8 j) turns into = Hvfct* ■ k. Hamil- 
tonian has eigenvalues: E' n k = n ■ fryp|k| = E n ^ that are degenerate with 
eigenvalues of Hamiltonian so that K' point represents only an additional 
degree of freedom like electron spin. In other words we can limit ourself to the 
behavior around K point if we note that each state is four fold degenerate i.e. 
two spin and two valley (K — K') degenerate. 



Finally, let us find the electron density and electron current density operators 
for Dirac electrons in graphene. To start, note that the electron momentum is 
p = Kk, which can be written as an operator in the coordinate representation 
p = — i/iV, so the Dirac Hamiltonian (I2.27P can be written as: H = —icr ■ V. 
If we now describe this Dirac electron by a wave function ip a (r) = (r\a), then 
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the electron particle density is simply p a = \ip a (r)\ 2 so the density operator is 

p op (r) = 5(t-t op ). (2.30) 

To find the electron current density we can apply the equation of continuity: 
— e §? + ^ ' J = (here we take e > so that — e denotes the electron charge), 
with an equation of motion Hip = ih^. This yields the electron current 
density: j Q = — euF^ a (r)*cr'0 Q ,(r), i.e. the current density operator: 

jop(r) = -ev F crS(r - r op ). (2.31) 

At last, the Fourier transforms of these quantities are given by 

Po P (q) = -^e-^ r -,and (2.32) 

jop(q) = -jZ<re**"*. (2.33) 



2.1.3 Electron-phonon interaction 

Since graphene is a 2D crystal with two atoms per basis, there are also two 
optical phonon branches (transverse and longitudinal) that are degenerate at 
energy ftwo = 0.196 eV and mostly independent of wave vector q (for long wave 
modes q « 2ix /a). Let us denote by u(R) = [u^(R) — ub(R)]/\/2 motion of 
the basis atom A relative to the atom B in the unit cell at the position R (see 
figure 12. 3p . If A and B where oppositely charged ions like in polar crystals, 
then their motion would result in the electric dipole moment i.e. electric 
field in the direction of the vector u and strong electron-phonon interaction. 
However, since A and B are completely equivalent carbon atoms, graphene 
belongs to the class of covalent crystals, and electron-phonon interaction is 
considerably reduced compared to the case of polar crystals. We will also 
see that electron-phonon interaction in graphene acquires unusual form in the 
vicinity of the Dirac (K and K') point and we will demonstrate that optical 
phonon oscillation creates effective electric field that is perpendicular to the 
vector u. That fact will lead to peculiar mixing of plasmon and optical phonon 
polarizations. 

The rigorous calculation of electron-phonon interaction in graphene is given 
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(b) 

f-K 




Figure 2.3: a) Atom motion in graphene during the longitudinal optical 
phonon oscillations, b) Motion of basis atoms described by a vector u in the 
real space induces vector potential A that moves Dirac points in the recip- 
rocal space. Dirac point symmetry causes unusual polarization of this vector 
potential: A _L u. 



in references pTl 153"] , while we only sketch here the main steps. Let us start 
with the tight binding Hamiltonian (I2.18P 



T(k) 
T*(k) 



(2.34) 



The effect of phonon on the electron motion can be simply found by consid- 
ering the change in the hopping integral (—70) with the change in the nearest 
neighbor distance. Let us now observe atom A at a position Ra and a neigh- 
boring atom B at a position Rb = Ra ~ t~i whose equilibrium relative distance 
is simply |t;| = b. If we move these two atoms out of equilibrium positions, 
new distance is: \ti + ua(Ra) — u b (Ra — Tj)|, and the leading order change 
in the hopping integral is 



-7 



-7o - 
-7o - 



db [l 
57o(6)l 
db b 



Ti + u A (R A ) - u B (R A - T t ) \ - b] 
n ■ [u a (Ra) - ub(Ra - Ti)] . 



(2.35) 



Since we are interested in long wavelength optical phonons (q << 2ii/a), 
then instead of discrete vector R, we can write a continuous coordinate r in 
the expression u(R) = [u^(R) — u^(R)]/\/2, i.e. we write 



ua(Ra) - u b (Ra - Ti) w u A (r) - u B (r - n) w u(r)\/2. 



(2.36) 



Finally the change in the hopping integral f l2.35p . in the long wavelength limit, 
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is given by: 

-7»-7o-^fr,.u(r). (2.37) 

Note here that all three neighboring carbon atoms (/ = 1,2,3) see the same 
phonon amplitude u(r), which will not be true in the case of finite wavevector 
q. However this change in amplitude will come with an extra factor q ■ a, so 
unless we are working with phonon wavevectors on the order of Brillouin zone, 
the long wavelength limit is a great approximation concerning the interaction 
between electrons and optical phonons. 



We can write phonon motion u(r) as a sum over normal modes 

u w = E 7/m Q ^ elCi ' r - (2 ' 38) 

Here M is the mass of a carbon atom, fj, — L,T denotes longitudinal i.e. 
transverse polarization, and if we define an angle <f(q) = tan -1 (q y /q x ), then 
polarization vectors are given by 

e qL = ^(cos f(q)x + sin f(q)y), and (2.39) 

e qT = i(— sin (p(q)x + cos (p(q)y). (2.40) 

Finally we can write the phonon amplitude Q qAt through the creation (fr q/ J and 
annihilation operators (b nfl ) as 

Q^ = ^(b m + bl m ). (2.41) 

To find how the phonon motion u(r) influences the electrons around the K 
point let us change the origin of wave vector as before: k — > k + K. Now the 
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function T(k) becomes 



T(k, r) = J2(-"?) Ze ~ i{k+K) ' Tl 
i 

- E (-0* - T T ' ■ u «) Ze " <<k+K,T ' 

- E - nlr x T < ' u(r ») Ze ~"' r, < 1 - !k ' T ' ] - (2 - 42) 

By looking at the leading order expansion in the phonon motion u(r) and 
electron wave vector k we have 



r( k ,r)^ 7oik -^fu(r)).^r,Ze— • 



(2.43) 



We recognize the first part of the expression (I2.43P from the equation ( I2.24p 
for bare Dirac electrons 

T (k) = 7o *k -J2-riZe- iKT < = ^-a lo {k x - iky), (2.44) 
i 

while the other part of the sum (12.431) gives 

T e _ pft (r) = — -u(r)-^T,Ze 

i 

^70(6) 72^3 . 

= ^ ^— a(-zn x -n y ). (2.45) 

With the substitution: hvp = 2 a ^°> expressions above transform into a sim- 
pler form 

T (k) = hvp{k x — iky), and (2.46) 

T e - ph {r) = Hv F —^-^(iu x + u y ). (2.47) 

Finally since T(k, r) = T (k) + ^.^(r), we can also write for the total Hamil- 
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tonian = + H e ~ph where 



o 




T (k) 
T*(k) 

T e _ ph \ _ fry (b) \/2 / Uy + iUz 

T e ^ ph i ~ VF db 6 7 o I u y - iu x 



(2.49) 



If we introduce here the notation cr x u = o^ity — o^Wx , then we can write 
equation (I2.49P in a convenient form [5TJ E2J 153] : 

# e _ p , = -fe^^ x u(r). (2.50) 
do 670 

From this expression we can immediately see the unusual property of the 
electron-phonon interaction in the vicinity of Dirac point. Namely the to- 
tal Hamiltonian in the presence of the phonons, can be obtained from the 
bare Hamiltonian = hvpcr • k, by a simple substitution: 

k x -> k x + Ku y , (2.51) 

k y -> k y - Ku x , (2.52) 

where K = ^jjj^ &^ • But this is precisely equivalent to the action of the vector 
potential A: 

hk ->■ hk + eA. (2.53) 

In other words influence of phonons on the electron motion is equivalent to 
the presence of vector potential with components: A x oc u y and A y oc — u x . 
This will in turn lead to the unusual mixing of plasmon and optical phonon 
polarizations. To understand this let us assume that the phonon wave vector 
is oriented in the y direction (q = qy) and let us look at the longitudinal 
optical phonon motion that have u x = and u y 7^ (see figure |2T3|) . Then 
the phonon influence is given by transverse vector potential since A x 7^ and 
A y = 0. In other words longitudinal phonon oscillation is equivalent to the 
transverse vector potential oscillation i.e. transverse electric field. On the other 
hand, since plasmons are collective charge density oscillations, accompanied by 
a longitudinal electric field, there won't be any interaction between plasmon 
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and longitudinal optical phonon. We will further show that there is a strong 
interaction of plasmon and transverse optical phonon which is a very counter 
intuitive result from the perspective of polar crystals. 

The simplest way to analyze the electron-phonon interaction in graphene is to 
show how the phonon amplitude couples to the electron current density. In that 
regards let us take electron-phonon Hamiltonian (12.501) and write expansion 
of the phonon motion u(r) over the normal modes from equation (12.381) to 
obtain: 

Hc _ ph = _ ft „ F ^) g_L= £ e> , V x e wQw . (2.54) 



Here we recognize the current density operator y = — ^f 1 cre* q ' r from equation 
( I2.33p . and if we introduce the factor F = f 6^ ^nm ' we can nnan V wr ite 
for the electron-phonon interaction Hamiltonian: 

H e _ ph = L 2 Fj2ii x e q/1 Q qM . (2.55) 

A more convenient way to write electron-phonon interaction is to show how 
the phonon amplitude couples to the electron density. In that respect, let us 
define the quantities: 

E^l = e q L • y + ie^L ■ % = i sin <p(q) — cos </?(q) = — e - ^^, and (2.56) 

EqT = eqT ■ V + ^qT ■ x = i cos y?(q) + sin <p(q) = ie~ lv> ^ . (2-57) 
Then by using the normal mode expansion (12.381) we obtain 

u y + iu x = , Qq/j(e qAt ■ y + ieqn ' x)e q = Q^E^e q . 

q,M q,M 

(2.58) 

Finally the electron-phonon interaction Hamiltonian ( 12. 49ft can be written as: 



If we now define 

<9 7o (fr) y/2 1 

g = hv F — - — ^=,and (2.60) 

db b-foVNM 
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qfj, — 




(2.61) 



then we can write electron-phonon interaction as a coupling between phonon 
amplitude Q m and electron density operator pt from equation (12.321) as: 



Formula (12.551) and (12.621) are equivalent. However, the response of the sys- 
tem to the interaction Hamiltonian (I2.55P is most easily described by utilizing 
the current- current response function while the response to the interaction 
Hamiltonian (12.621) is most easily described by the density-density response 
function. 

2.2 Response functions 

In this section we present response functions which describe response of our 
system to an external perturbation. Specifically we calculate graphene's con- 
ductivity, density- density, and current- current response functions in the weak 
coupling approximation i.e. linear response theory. These three functions are 
all connected by simple relations, however it will be more convenient to use 
one or another depending on the specific nature of the problem being studied. 

2.2.1 Conductivity 

Semiclassical model - Drude conductivity 

If we are only interested in the response of the graphene under the influence 
of external electromagnetic field, we can simply calculate the conductivity 
function a(u). The semiclassical model gives a simple relation for the Drude 
conductivity [38] 





(2.62) 





(2.63) 
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spin and two valley degeneracy. Semiclassical model is simply a generalization 
of the Drude model for free electrons to the case of an arbitrary band structure 
E(k), however we will see that it can describe lot of interesting phenomena in 
a qualitatively correct way. At zero temperature one has —df/dE = 5(E — fi) 
and it is straightforward to show that for the case of Dirac electrons in graphene 
H k = hvpcr ■ k, the Drude conductivity is given by 

7TAT 1/r — IUJ 

It is a slightly more tedious task to show that at finite temperature T one has 

a D (oj) = e 2 —\n 2 cosh -£-) -. 2.65 

7r/r V 2A;T/ 1/r — ioj 

Fermi's golden rule - interband conductivity 

The semiclassical model has a serious limitation since it cannot describe transi- 
tions between different bands which is particularly important in graphene 
that has zero band-gap between valence and conduction bands. To take into 
account these interband transitions we will calculate the response of graphene 
to an external electric field, in the first order perturbation theory (using the 
Fermi's golden rule). 

Let us imagine that an electromagnetic plane wave of frequency u is incident 
under the normal angle onto the graphene sheet. We can choose the gauge so 
that the scalar potential <p = 0, while the vector potential A = Aoe~ tu)t , 
so that the electric field is given by E = — ^ = iuA = E e _JaJ * , and 
Eo = ioj Aq. Electrons in graphene are described by a Dirac Hamiltonian 
( I2.27P Hq = Hvf<t ■ k = vf<t ■ P, where p is the graphene's electron momentum 
so that the interaction with the vector potential is simply described by a sub- 
stitution p — > p + eA. In other words, the total Hamiltonian in the presence 
of electromagnetic field can be written as H = vp& ■ (p + eA) — Hq + Hi nt , 
where the interaction part of the Hamiltonian is given by 

H in t = ev F o- ■ A = —a- ■ E e~ iujt . (2.66) 

Here we have kept only the time dependent part (e~ tuJt ) responsible for the 
absorption process. Then the Fermi's golden rule [BB] gives the probability for 
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(2.64) 



a transition from an initial state i to the final state /, with an absorption of a 
photon: 

rlin- . p 9tt 

(i\H int \f)\ 2 5(huj if - hw)f t (l - //)■ (2-67) 



dt h 

The total power absorbed from the incident wave can be written in two ways. 
First, one can write 

P ° = £^f- ( 2 - 68 ) 

On the other hand, since = er(u;)E(u;), one can write (for harmonic fields) 



P a = 25ft J j(w) • E*(oo)dr = 2Ka(w)|Eo| 2 L 2 , (2.69) 

where L 2 is the area of graphene sheet, and we used the fact that E(cj) = Eo is 
uniform along the graphene plane for the case of normal incident wave. Finally 
we have 

= 2^T7^ £ ^\{i\H int \f)\ 2 S(hw if - fu^Ml - f f ). (2.70) 

Now, let us denote the initial (final) state of the electron by a band index n 
{p!) and a wave vector k (k') i.e. \i) = \nk) (\f) = |n'k')). Without loss 
of generality we can assume that the electric field is polarized along the x 
direction: E = Eq'x. Then we can write for the matrix element: 

(i\H int \f) = —E (n'k'\a x \nk}. (2.71) 

Further on, by using explicit form (I2.29P for the Dirac electron wave function 
^> n ,k, it is simple to show that 

(n'k»k> = J < /)k ,(r)^n,k(r)dr = ^^(ne"^ + rie**), (2.72) 

so we obtain expressions for the matrix element 

|(*|# mt |/>| 2 = ^|£ | 2 5 klk i(l +rm' cos 20 k ), (2.73) 
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and conductivity 

tor(u) = ^4Vl / kdk f d6^^ihl + nn'cos26 k )x 
2 ^^4ir z J I n ur 2 

n,n' 

5(n'hv F k - nhv F k - hu)f nk (l - f n > k ), (2.74) 

where we also took into account 2 spin and 2 valley degeneracy. We now 
take into account only (interband) transitions between conduction and valence 
bands, because the intraband transitions are already taken into account by the 
Drude conductivity. After lengthy but straightforward calculation, one obtains 
simple expression for the real part of the conductivity 

= ^/(-«w/2)[l - /(fiw/2)]. (2-75) 

It is instructive to look at this result at zero temperature 

e 2 

ftafa) = ^0(hj - 2/i). (2.76) 

Here 6(x) is a simple step function [8(x < 0) = and 6(x > 0) = 1]. The real 
part of the conductivity provides us with absorption of the electromagnetic field 
incident on a graphene sheet. We see that there is no absorption for faj < 2/i 
which is result of the Pauli exclusion principle. On the other hand above this 
threshold, when ftjjj > 2/i one will have uniform absorption. Since the incident 
energy flux is given by Wi = 2\E \ 2 /fi c (see reference [SZj)j and the absorbed 
energy per unit time per unit area is given by W a = P a /L 2 = 2^R,a(uj)\Eo\ 2 (see 
equation (I2.69P ). the absorption coefficient can be written as 

2 _ W a _ /i ce 2 _ 
\a\ - — - — -2.3%. (2.77) 

This result has been confirmed by experiment [57] • Further on, note that if we 
include the emission process, then we obtain the following expression for the 
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conductivity: 



e 2 



tor(u) = -[f(-hoj/2)[l - f(fkj/2)] - [1 - f(-hw/2)]f(hw/2)] 

= ^\f(-fa/2)-f(tiu>/2)]. (2.78) 

Finally we can obtain the imaginary part of the conductivity by using the 
Kramers- Kronig relations [56] : 



2u; f°° 9far(u/) 



= V / -fr^duf 



e 2 4fkJ f°° f(-e) - fit) , 

-v I ; \ 2 /fc\o rfe - ( 2 - 79 ) 



4fr 7r 7 (2e) 2 -(^;) 2 
It is convenient to introduce the following function: 

KT + CO&ii fcT 

Then, we can simply write for the total interband conductivity (see also [61]): 



<"^ = Th{ g ^ + '—1 ww*)' (2 ' 81) 

In the last expression, we took into account that principal value of the integral 
with G(hw/2) equals to zero, which removes singularities from the integral in 
the imaginary part of the conductivity. 



2.2.2 Density-density response function 

We now proceed to a more formal, but powerful, aspect of linear response 
theory by looking into the density- density response function. In the last section 
we assumed that there is no spatial dependence of external perturbation and 
calculated only frequency dependence of the conductivity. Let us assume that 
graphene is placed in an external scalar potential of arbitrary spatial and time 
dependence 

cp ext (r,t) = I e-^duYs^VexM^). (2.82) 
J q 
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Now, scalar potential simply couples to the electron charge density so one can 
write the interaction Hamiltonian [SB] 

H mt = I e-^duL 2 J2(-^e X t(<l^))pl (2.83) 
J q 

We now assume the weak coupling between the system (electron density) and 
a probe (external potential) so that we can focus on a single (q, u>) component. 
The induced electron particle density is then given by 

(/w(q,o;)) =x(q.w)(-ev7« BB t(q,w)), (2-84) 

where the density-density response function is given by [56] 

^ P -P E b / 1 1 \ 

X (q, uj)=L^ (a|pi|6> I 2 t ^—T - " ^ Z1T + • 

Z H \nuj — njjj ab + iT) rujj + nw ah + iri I 

a.b ' 

(2.85) 

Here Z = e~^ Eh is the partition function, and hw a b = E a — E^. Further on 
| a), and E a are exact many body state, and energy of the system in the presence 
of the perturbation. In other words we can write H\a) = E a \a) where H = 
Hq + Hi n t is the total system Hamiltonian given by the sum of the Hamiltonian 
in the absence of perturbation (Hq) and the interaction term (H int ). Equation 
02. 85!) is exact in the limit of weak coupling (i.e. linear response), however one 
first needs to find the exact eigenstates of the total Hamiltonian H which is 
not an easy task. We shall deal with this issue by working in the self-consistent 
approximation i.e. by introducing simple, yet powerful, concept of screening. 
In that regard let us note that the induced charge density (pmd(<l, w)) will 
be accompanied by an scalar potential <£w(q, w) which can act back on the 
electrons through the interaction Hamiltonian f)2.83p . In other words, instead 
of equation (12.84ft we should write the self-consistent equation for the total 
induced particle density 

(p ind (q, u)) = x(q, w)(-ey? exi (q, u) - e(p ind (q, cu)). (2.86) 

However, x(q, u) is now the screened density- density response function which 
is again given by the equation (12.851) . only \a, b) are now simply the eigenstates 
of the noninteracting Hamiltonian H . This is the lowest order approximation 
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which can also be traced down to the random phase approximation. For a 
system of Dirac electrons, described by a wave functions VVi,k given by equation 
(I2.29p . one then obtains for the screened response function |56j : 

X(q, -) = £ I <-' k + q|e^>k)| 2 - f ^ 

Li KM — -^n'k+q + -&nk + ^ 

nn'k 

nn'k 



We will be particularly interested in the dielectric function of this system so 
we need to find the relation between the scalar potential <fi n d and the induced 
surface charge density —e(pi n d(r,t)) = — e(pi n d(q, u))e l<i ' r e~ luJt . Let us define 
here the vector r = xx + yy which lies in the graphene plane (located at 
z — 0) while z axis is perpendicular to graphene plane. Further on we assume 
graphene is sitting in between two dielectrics of permittivities e r \ (z < 0) and 
e r2 (z > 0). If we work in the electrostatic approximation (q » u>/c) then 
the scalar potential induced by the surface charge density located at the plane 
z = is simply given by 

<p ind (r, z, t) = <^ nd (q, uy^-^e- 1 ^. (2.88) 

The electric field is given by E = — V</>. We can now separate the electric field 
E = E r + E z z into component along the graphene plane E r = — V r <y9 which is 
given by expression 

E^V, z, t) = -tp ind (q, c^qe^-^e-^, (2.89) 

and component perpendicular to the graphene plane E z = —dip/dz which is 
given by expressions 

Ef*(r, z>0,t) = ^ nd (q, u)iqe^ z , (2.90) 

ErV, z<0,t) = -<p ind (<l, u^qe^^e-^. (2.91) 
Further on, the Gauss law can be written as a boundary condition across the 
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graphene plane as [67] 



-e(p ind (r, t)) = [B md (r, z = 0+ t) - B md (r, z = 0- t)} ■ z 

= e e rl E™ d (r,z = + ,t) - e e r2 E™ d (r,z = 0" i). (2.92) 

Then by using the decomposition into Fourier components and equations ( I2.90p 
and f)2.9ip we obtain desired relation between the induced charge density and 
corresponding induced scalar potential: 

- e(/w(q, u)) = q(pind(q, oj)2e r e . (2.93) 

Here e r = (e r i + e r i)/2, and we can introduce the external charge density 
corresponding to the external potential by the same relation 

- ep ext (q, u) = qtpext(<l, uj)2e r e . (2.94) 

Let us now define the graphene dielectric function e(q, u) as [56] : 

e(q,co) p ext (q,u) 



e r Pext(q, of) + (pmd(q, u)) ' 
Then from equations (I2TH5]) . (I2T9"BD and (I2~M1) we obtain 



(2.95) 



^=l-5^-x(q,«). (2.96) 

Note that the zero of dielectric function (e(q, u) = 0) defines the collective 
electron oscillation (plasmon) which is the core subject of this thesis. 



Finally let us find the relation between density-density response function 
x(q, w) and conductivity er(q, ou). If we introduce the total scalar poten- 
tial iptot = fext + find, then by using equation ( 12. 86ft . we can write the in- 
duced surface charge density as — e(p in d(q, u>)) = x(q, u)e 2 ip tot (q, u>). On the 
other hand Ohm's law gives the induced surface current density (jind(q, oo)) = 
er(q, w)E*°'(q, u) , while the electric field can be found from equation (12.891 ): 
E*°*(q, cu) = —cp tot (q,cu)iq. Finally, equation of continuity can be written with 
Fourier components as — e(p in d(q, uj)) = q- (jind(q> so we obtain desired 
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relation: 

a(q,u) = i— x(q,w). (2.97) 

Note however that cr(q, ui) refers only to the longitudinal conductivity since 
the scalar potential alone is not enough to decribe the transverse fields. 



2.2.3 Current-current response function 

In the last section we described response to the external scalar potential which 
we now supplement by calculating response to the external vector potential. 
Let us then start with the Hamiltonian ( 12.271) describing free Dirac particles: 
Hq = vf<t ■ p, where p is the electron momentum. In the presence of external 
vector potential A ext (r,t), one can write for the total Hamiltonian H = vpcr ■ 
(p + eA ext (r, t)) = H Q + H int , where the interaction part of the Hamiltonian 
is given by: H int = evpcr ■ A ext (r,t). We can now decompose vector potential 
into Fourier components to obtain: 

H int = £ e^ r ev F( r ■ A ext (q, t), (2.98) 
q 

then by using the current density operator from equation ( I2.33|) we can write 

H int = -L 2 ■ A -t(q, <)■ (2-") 
q 

It is now convenient to introduce the longitudinal (Vl = V-e* L ) and transverse 
(y T = V ■ e* T ) vector components by using the polarization vectors from 
equations f!2.39j) and (I2.40p . We can now write the interaction Hamiltonian 

H int = -L 2 ^^ ■ A ext ^(q, t) 
q,M 

= J e-^duj(-L 2 ) ■ ^Wq,"). ( 2 - 100 ) 

q^ 

Finally, by assuming the weak coupling between the external probe and our 
system, precisely like in the last section, we obtain the induced current density: 

(jindA^u)) = x^^i-AxtA^^))- (2.101) 
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Here the current- current response function is given by [56] 

^ p-PEb / 1 1 \ 

S—fZ \hW — tiWab + 11] hw + nWafe + ^ / 

(2.102) 

We can now use the free electron states to write the screened function: 

X ,(q,a,) = |(n'k + q|jt |nk>| 2 - ^ " ^ - ■ (2.103) 

nn'k 

At last, by using the exact form of the electron wave function from equa- 
tion f)2.29p . we obtain different expressions for the longitudinal and transverse 
current-current response functions: 

Xl(q,w) = — 2-4^-[l+nn cos(# k + # k+q )]— 

n^k ~~ -^n'k+q + Alk + 

(2.104) 

XT(q,u;) = ^4 £ i[l - nn' cos(0 k + k+q )]- ^ k " ^ - . ■ 

nn'k 

(2.105) 

Note here that expressions (I2.104p and (12. 104j) actually diverge if we use 
Dirac states ( I2.29P instead of actual electron states in graphene limited by 
some band cut-off. However, his subtlety can be easily solved by subtracting 
from Xi(q,w) [xt(<1,u)] the value Xl(<1,w = 0) [%y(q ->• 0,u = 0)] to take 
into account that there is no current response to the longitudinal [transverse] 
time [time and space] independent vector potential, see [511 ES] for details. 

Let us also find relation between the conductivity and the current-current 
response function. Note that the electric field is given by E = —dA/dt so that 
E(q, uS) = zu;A(q, cj). Then we can write equation ( I2.10ip as (j ind ^(q,u)) = 
Xu(q, u)^E exttfl (q, oj). In other words desired relation is simply: 

i 

cr u (q,w) = -Xn(q,w)- (2.106) 

Note here that longitudinal conductivity 0z,(q, ui) (describing response of a 
system to the longitudinal field) is generally different from the transverse con- 
ductivity OY(q, u>) (describing response of a system to the transverse field), 
unless we are working in the limit of small wave vectors (q — > 0). 
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2.2.4 Fluctuation-dissipation theorem 

In this section we derive relation between current- current correlation function 
and the current-current response function at finite temperature, which is given 
by the fluctuation-dissipation theorem. We will use this result later to calculate 
the radiative heat transfer between two graphene sheets. 
We start with the current- current correlation function: 

K,(r,t,v',t') = (j,(v,t)jl t (r',t')). (2.107) 

Due to translational invariance in space and time we can write K^r, t, r', t') = 
K^r — r',t — t') = K^(d,T) , where we have denoted by: d = r — r' and 
r = t — t'. Then the Fourier transforms from the space and time domains are 
respectively given by 

K^t) = ^J ^(d,r)e-^ d rfd, (2.108) 

K,(q,u) = ^ J K,(q,r)e WT dr. (2.109) 

It will be more convenient for us to use these relations in a slightly different 
form. In that regards let us use relations (I2.108P and (12. 109[) with translational 
invariance in space and time, respectively, to show that 

0,(q)4(q')> = 71 // (U^li^e-^e^dvdY 1 

= j2 J e^'-^'dr'-^ j K^r - r^e^^dr 

= <w^(q), (2-110) 



<J»J>')> = J {j,(t)jl(t'))e^e-^'dtdt' 

= — f e^-^'dt'— [ KJt - ty-^-^dt 

= 5(uj-u')K^(uj). (2.111) 

Relations ( I2.110p and (12.1 lip simply state that there is no correlation between 
different q or different u components. We can join these two relations in a 
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single one 

(i ( .(q I w)i!(q'y)) = W(w-^(q,w). (2.112) 

To find the K^(q,u) let us note that evolution of current operator, in the 
Heisenberg picture, is given by j M (q, r) = e tHT ' h jJ l (q,0)e~ iHT ' fi . Now we can 
write 

^(q,r) = (j M (q,r)jt(q,0)) 

= (e^ T/rt j'J(q,0)e"^jt(q,0)) 

= E -^( & l^ T/ %(q)l«)(«|e-^'; t (q)l&) 



0,6 



E i 7^i< a i4(q)i 6 >i ae * ,WT - ( 2 - 113 ) 



a, 6 



Finally, the Fourier transform of this expression is given by 

-j3E b 



^(q,^) = E^-|Hj:(q)l&)| 2 ^-^). (2.114) 



z 

a,b 

Note however that the imaginary part of the response function, calculated in 
equation (12.102]) . is given by 

^(q,o;) = L 2 E -^Mjt(c0\b}\ 2 — [5(co - u ab ) - 5(u + u ab )\. (2.115) 

a, b 

We immediately see that correlation function is related to a response function 
in a simple manner: 

3 XM (q, w) = -^L 2 [i^(q, w) - ^(q, -w)] . (2.116) 

By applying the detail balancing condition here, we can write 

%XM^) = -\L 2 [l-e-^] Kfau). (2.117) 

Finally we have 

K»(d,u) = —Y^^^XMM, (2-118) 
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or if we use the relation Xni^l^) — ~ wj<7 M (q, we can write this in a more 
convenient form as 



l hu i 



3fcr M (q,o;). 



(2.119) 



7T 1 - L 2 



This is in fact the well know fluctuation-dissipation theorem stating that the 
correlation function (K^) due to thermal fluctuations is directly related to 
the dissipation in the system (Sftcr^ or Sfrx^). This result will be of use in the 
following section. 

2.3 Radiative heat transfer 

In this section we analyze the radiative heat transfer between two graphene 
sheets separated by a distance D and held at temperatures 7\ and T 2 (see 
figure |2~4|) . To calculate the heat transfer we shall start by looking into cor- 
relations between electric currents induced by the thermal fluctuations in the 
first graphene sheet. Following that we shall use Green function technique to 
find the electromagnetic fields in the second graphene sheet, induced by the 
fluctuating currents from the first sheet. Finally heat transfer can be found 
by calculating Ohmic losses, induced by this electromagnetic field, within the 
second graphene sheet. 

In the last section we calculated current-current correlation function due to 
thermal fluctuations. Fluctuation-dissipation theorem (12.1 19j) and equation 
(12.1121) give the correlation function of the fluctuating currents in the first 
graphene sheet: 



To find the electromagnetic fields induced by these fluctuating currents we 
can use classical electrodynamics so we shall start with classical quantities and 
return to the quantum values only later when necessary. Since the system is 
translational invariant we can focus on a single q, w-component and write the 
Fourier transform of the surface current density from the first graphene sheet 
as 



(ji/,(q, ^)jj M (q,w')) = 5(uj- u') 



l hu i 



3&ri M (q,a/). 



(2.120) 



7r 1 - e-ft*" L 2 




q 



(2.121) 
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Further on, let us assume the most simple case where there is only vacuum in 
between and around graphene sheets. Then the electric field satisfies a simple 
wave equation 

(V 2 + oo 2 /c 2 )E(r,z) = 0. (2.122) 

This equation has a plane wave solution E(r, z) = E(q, w)e 4q ' r+nz , where we 
took into consideration that the periodicity in the xy direction is determined 
by the wave vector q. In other words we can write for the total wave vector: 
w = q + 7z, while the equation (12. 122ft requires: w 2 = |w| 2 = q 2 + 7 2 = u 2 /c 2 
i.e. the z component of the wave vector w is given by: 

7 = ^Vc 2 - q 2 . (2.123) 

Further on, since there is no free charge around graphene sheets, the Gauss 
law states that V • E(r, z) = 0. This means that E(q, uj) ■ w = i.e. electric 
field is transversely polarized so it is convenient to introduce unit vectors s 
and p that are perpendicular to wave vector w: 

s = q x z, and (2.124) 

p = w - 1 {--fq + qz). (2.125) 

In this way (s, w, p) is a set of right-handed orthonormal triad (see figure 12.4ft 
where 

w 

w = — = uT 1 (gq + 7z). (2.126) 

w 

We also note that there is a simple connection with the longitudinal and trans- 
verse wave vectors introduced before in this chapter: e q ^ = iq, e q ^ = —is. 

To match the boundary conditions given by the surface current density 
ji(q, oo) from the first graphene plane (z = —D) in the presence of the second 
graphene sheet (z — 0), we use the Green function technique from reference 
[60] which is particularly convenient for the layered structures like ours. In 
that manner one obtains different electric field component E(q, uj) depending 
whether we are located below the first graphene sheet (z < —D), in between 
the sheets (-D < z < 0), or above the second sheet (z > 0). Since we are 
interested in the field in the second sheet (2 = 0), it is easiest to look into the 
expression for the field above the second graphene sheet (z > 0) where one 
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(a) 



< s ggigsa 



7*Z% 



x5% 



-D-- 



(b) 




Figure 2.4: (a) Schematic diagram of the radiation transfer problem: a free 
standing sheet of graphene at temperature T\ is radiating to another free stand- 
ing graphene sheet at temperature T 2 and distance D away, (b) Polarization 
vectors defined in the text. 



obtains 



E(q, W ) = --^(sT^s + pT^p) •j 1 (q,u;). (2.127) 



Here we have explicitly separated s and p polarizations which have very differ- 
ent behavior, and T 12 is a transmission coefficient for a system of two parallel 
graphene sheets given by [60] 

tit 2 e^ D 

T i2 = i 2nD - 2 - 128 

1 — r\r 2 e 1 

Note that the same expression is valid for s and p polarization, but reflection 
r and transmission t coefficients are different for different polarizations. It is 
a simple manner of elementary electrodynamics to demonstrate that these are 



LOUT 

2-yeoc 2 
1 + 



1 I LOUy ' 

1 ~+~ 



7°X 



(2.129) 
(2.130) 



r* = -f§^,and (2.131] 



tP = T-^T- ( 2 - 132 ) 
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Note also that transverse conductivity (<tt) determines the s-polarization, and 
longitudinal conductivity (cr^) determines the p-polarization. Finally the total 
electric field in the space above the second graphene sheet (z > 0) is given by 

E(r,z,t)= I e~ iult duJ2e ici - r+iwz E(q,uj). (2.133) 
J q 

So the field, precisely at the second sheet [z = 0) is 

E(r, t) = E(r, z = + , t) = ! e~ iwt du ^ e iq ' r E(q, cu). (2.134) 

J q 



At last, the heat transfer from the first graphene sheet to the second graphene 
sheet is simply given by Ohmic losses induced by this electric field. The power 
dissipated per unit area is given by [67] 




dudu'e-^-^ J2(q, w) • E*(q, u'). (2.135) 



Let us now take into account that current densities ji,2(q> oS) have only vector 
components along the graphene (xy) plane. Then due to equation (I2.125P one 
has p • ji(q, uj) = — • ji(q, oj), and we can write equation (12.1271) again as 

E(q,c) = (sT^S+ (^q- ^ T f 2 q) ■ j^a,). (2.136) 

Further on, the scalar product in the equation (I2.135P can be written as j2-E* = 
j2 ■ E*, where E r = E — (E ■ z)z is the projection of the electric field vector to 
the graphene (xy) plane. In that way we can write equation (12.1361) as 

E r (q,w) = fe 2 s+ ^q^q) ■ ji(q,w). (2.137) 

2e c^7 \ w A J 

Let us note here again that q = ie* L and s = — ?e* T while the longitudinal 
(/i = L) and transverse (fi = T) components of the current density are defined 
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as: jn(tl,u) = j(q, u) ■ e*^. Then we can write equation (12. 137j) as 

Erfaw) = (-< T Tf 2 j 1T (q,a;) + -^e^Tf^q, u)j . (2.138) 

Finally due to Ohm's law j 2A1 (q, cu) = o" 2/i (q, u)E riM (q, oj) we have 



U) 2 

(j 2 (q,w) -E*(q,w')) = . 2 4I a 2T (q,u;)(j 1T (q,a;)j 1T (q,u; / )*)|yi2| 2 



h-4La 2L (q,a;)(j li (q,a;)i li (q,a; / )*)|rf 2 | 2 . (2.139) 



Here we have explicitly written the ensemble average which requires us to 
calculate precise quantum correlations of the current density operator. We 
have also used relation 02.112]) which states that there is no correlation between 
different u components due to translational invariance in the time domain. In 
fact the current-current correlation function (12.1201) is given by 

(ji„(q, tt)jUq, = *(" - i _ ^ jfb^q, u). (2.140) 

Since the final result has to be a real quantity, we can simply look into real 
part of the expression (12. 139p 

• KM) - ^') l -T^^h x 

( 4e2 ^ 4 , |2 ^2T(q, co)^a 1T (q, u>) \T° 2 \ 2 + 



+ jT-2^2L(q,u)^r 1L (ti,u)\Tf 2 \ 2 ). (2.141) 
4e a; 



This can be written in a more transparent form by using reflection and trans- 
mission coefficients (12. 129ft - (12.132ft . However, since 7 = \/oo 2 /c 2 — q 2 we have 
to distinguish between the case of propagating waves in the far field (u/c> q) 
and evanescent waves in the near field (u/c < q). In the first case (co/c > q) 
one has 

— %—^a T = ' J ' , and (2.142) 

2e c 2 | 7 | 2\t s \ 2 

JtI 

2e uj 



1 - 


r p 


2 






2 


2\tP 


2 



^ = 0L2 ■ ( 2 - 143 ) 
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It is convenient here to define the following quantities 



oJ 2 



h s ff ((l,u) = 4e 2 c 4| 7 | 2 ^ aiT ^' ^)3frx2r(q, u)\T£ 2 \ 



(1- 




2 -|^i| 2 )(l- 
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I 2 - 1*2 
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and (2.144) 
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|2 



4e o; 



(1- 
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2 -l^l 2 )(l- 
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I 2 - 1*5 
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4|1 


'l'2 e 1 


2 



(2.145) 



where we have used expression (I2.128P for the transmission coefficient T\i- In 
the second case (cu/c < q) one has 



^ = w (2 ' 147) 



KMM = 4e 2 c 4| 7 | 2 ^(q^)^T(q,u;)|r i s 2 | 2 = p— — — z^p and 

(2.148) 

U|2 CV T .PCV r P„-2| 7 |D 

*W q ' w) " 4^^(q^)^(q,a;)|Tf 2 | 2 = }1 J^.^ - (2.149) 

At last we obtain for the heat transfer (equation (12.135]) ) from the first graphene 
sheet to the second graphene sheet 

ifl->2 = Hl-+2Jf + Hi-+2,nf> (2.150) 

where the far field (u/c > q) and near field (u/c < q) contributions are 
respectively given by 

H^jf = - I du - _ ^ pihw — tffiq, a), and (2.151) 
= \\ ^YZ^^KM")- ( 2 - 152 ) 
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In the same manner one can calculate heat transfer from the second graphene 
sheet to the first graphene sheet H 2 ^i, so the total heat transfer (H = Hi^ 2 — 
i?2->-i) between two graphene sheets can be written as 

H = H ff + H nf , (2.153) 

2 roc poj/c 

#// = -/ MQiu^-e^Tz)}—— / 2KqdqY^h» f {q,u), (2.154) 
71 Jo K 171 ) Jo 

cy pOO i pOO 

H nf =- cL;[e(o;,T 2 )-e(a;,r 1 )] 7 — - / 27rgd g yV / (<z, w). (2.155) 

^ •/ l Z7r J ./w/c 



Here we have introduced the Boltzman factor: Q(u,T) = hid/ [er^ — 1), which 
comes about since the zero point energy cancels when taking the difference 
between emission and absorption. We write here again functions and 
for the sake of clearance 



(1- 




2 




' 1 


2 )(1- 
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4|1 




2 



h fM> ^ s l \u_J^Dv > ( 2 - 156 ) 



CV^C^./* -2| 7 |£> 

KM^ = {l _\»X-^ D r (2 ' 157) 

Note that for the case of black body which has perfect absorption 1 = \a\ 2 = 
1 — |r| 2 — \t\ 2 , i.e. zero reflection or transmission (r — t — 0), equation (I2.154p 
simply gives the Stefan-Boltzman law: 



To summarize, in this section we have calculated the total heat transfer, 
that is, the transfer of heat energy per unit time per unit area between two 
graphene sheets at different temperatures. Total heat transfer H = Hff + 
H n f has a contribution from the propagating waves in the far field {Hff) 
and evanescent waves in the near field (H n f), given by equations (I2.154p and 
(I2.155p . respectively. 
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Chapter 3 



Plasmonics in graphene 



In this chapter we investigate plasmons in doped graphene and demonstrate 
that they simultaneously enable low-losses and significant wave localization 
for frequencies of the light smaller than the optical phonon frequency hjJoph ~ 
0.2 eV. Interband losses via emission of electron-hole pairs (1 st order process) 
can be blocked by sufficiently increasing the doping level, which pushes the 
interband threshold frequency Co>j nter toward higher values (already experimen- 
tally achieved doping levels can push it even up to near infrared frequen- 
cies). The plasmon decay channel via emission of an optical phonon together 
with an electron-hole pair (2 nd order process) is inactive for uj < uo P h (due 
to energy conservation), however, for frequencies larger than uo P h this decay 
channel is non- negligible. This is particularly important for large enough dop- 
ing values when the interband threshold cjj nter is above tooph- i n the interval 
ujoph < u) < Umter the 1 st order process is suppressed, but the phonon decay 
channel is open. In this chapter, the calculation of losses is performed within 
the framework of a random-phase approximation (RPA) and number conserv- 
ing relaxation-time approximation [31]; the measured DC relaxation-time from 
Ref. [5] serves as an input parameter characterizing collisions with impurities, 
whereas the optical phonon relaxation times are estimated from the influence 
of the electron-phonon coupling [35] on the optical conductivity [36J. 

In Sec. 13.11 we provide a brief review of conventional surface plasmons and 
their relevance for nanophotonics. In Sec. 13.21 we discuss the trade off between 
plasmon losses and wave localization in doped graphene, as well as the optical 
properties of these plasmons. We conclude and provide an outlook in Sec. 13.31 
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3.1 Surface plasmons 



frequency [THz] 
600 550 500 450 

--- Re[q]/Im[q] / j c j 
10 X ' k i X *e 1/ 

8.5 0.55 0.6 0.65 

Figure 3.1: (a) Schematic description of a surface plasmon (SP) on metal- 
dielectric interface, (b) SP dispersion curve (solid blue line) for Ag-Si inter- 
faces; dotted blue is the light line in Si; dashed red line denotes the SP reso- 
nance, (c) Wave localization and propagation length for SPs at Ag-Si interface 
(experimental Ag losses are taken into account). 




Surface plasmons (SPs) are electromagnetic (EM) waves that propagate 
along the boundary surface of a metal and a dielectric [see Fig. 13.1( a)]; these 
are transverse magnetic (TM) modes accompanied by collective oscillations of 
surface charges, which decay exponentially in the transverse directions (see, 
e.g., Refs. [T2~l ITB] and Refs. therein). Their dispersion curve is given by: 

u I e r e(u) 

Qap=-\ — —r^ 3- 1 

[see Fig. 13.1T b)]: note that close to the SP resonance (u = ujsp), the SP wave 
vector [solid blue line in Fig. 13.1( b)] is much larger than the wave vector of the 
same frequency excitation in the bulk dielectric [dotted blue line in Fig. 13.1( b)], 
As a result, a localized SP wave packet can be much smaller than a same 
frequency wave packet in a dielectric. Moreover, this shrinkage is accompanied 
by a large transverse localization of the plasmonic modes. These features 
are considered very promising for enabling nano-photonics [121 H31 HH 1X5] . 
as well as high field localization and enhancement. A necessary condition 
for the existence of SPs is e(cu) < —e r (i.e., e{u) is negative), which is why 
metals are usually used. However, SPs in metals are known to have small 
propagation lengths, which are conveniently quantified (in terms of the SP 
wavelength) with the ratio ^R,q sp /^sq sp ; this quantity is a measure of how many 
SP wavelengths can an SP propagate before it loses most of its energy. The 
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wave localization (or wave "shrinkage") is quantified as \ a ir/X sp , where X a i r = 
2hc/uj (the wavelength in air). These quantities are plotted in Fig. 13.1( c) for 
the case of Ag-Si interface, by using experimental data (see [H] and references 
therein) to model silver (metal with the lowest losses for the frequencies of 
interest). Near the SP resonance, wave localization reaches its peak; however, 
losses are very high there resulting in a small propagation length I ps 0.1 X sp ~ 
5nm. At higher wavelengths one can achieve low losses but at the expense of 
poor wave localization. 

3.2 Plasmons and their losses in doped graphene 

Graphene behaves as an essentially 2D electronic system. In the absence of 
doping, conduction and valence bands meet at a point (called Dirac point) 
which is also the position of the Fermi energy. The band structure, calculated 
in the tight binding approximation is shown in Fig. 2(b); for low energies the 
dispersion around the Dirac point can be expressed as E n ^ = nvFh\k.\, where 
the Fermi velocity is vp = 10 6 m/s, n — 1 for conduction, and n = —1 for the 
valence band. Recent experiments [37] have shown that this linear dispersion 
relation is still valid even up to the energies (frequencies) of visible light, which 
includes the regime we are interested in. 

Here we consider TM modes in geometry depicted in figure [3T21 (a), where 
graphene is surrounded with dielectrics of constants e r \ and e r 2- Throughout 
the paper, for definiteness we use e r \ = 4 corresponding to SiC>2 substrate, and 
e r 2 = 1 for air on top of graphene, which corresponds to a typical experimental 
setup. TM modes are found by assuming that the electric field has the form 



After inserting this ansatz into Maxwells equations and matching the boundary 
conditions [which include the conductance of the 2D graphene layer, a(u,q)}, 
we obtain the dispersion relation for TM modes: 



E z = Ae iqz ~ QlX , E y = 0,E X = Be iqz ~ QlX , for x > 0, 
E z = Ce iqz+Q2X , E y = 0, E x = De iqz+Q2X , for x < 0. 



(3.2) 




(3.3) 
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Figure 3.2: (a) Schematic of the graphene system and TM plasmon modes. 
Note that the profile of the fields looks the same as the fields of an SP [Fig. 
13.1( a)]. (b) Electronic band structure of graphene; to indicate the vertical 
scale we show the Fermi energy level for the case Ep = 1 eV. (c) Sketch of the 
intraband (green arrows) and interband (red arrows) single particle excitations 
that can lead to large losses; these losses can be avoided by implementing a suf- 
ficiently high doping, (d) Plasmon RPA and semiclassical dispersion curves. 
Black solid (RPA) and black dot-dashed (semiclassical) lines correspond to 
£ri = £r2 = 1; Blue dashed (RPA) and blue dotted (semiclassical) lines corre- 
spond to e r i = 4 and e r 2 = 1. The green (lower) and rose (upper) shaded areas 
represent regimes of intraband and interband excitations, respectively. 



By explicitly writing the dependence of the conductivity on the wave vector 
q we allow for the possibility of nonlocal effects, where the mean free path of 
electrons can be smaller than q~ x [38]. Throughout this work we consider the 
nonretarded regime (q ^> uj/c), so equation (13. 3 p simplifies to 

q^Qi^Q 2 ~ e = - f v 3.4 

2 <T[u,q) 

Note that a small wavelength (large q) leads to a high transversal localization 
of the modes, which are also accompanied by a collective surface charge os- 
cillation, similar to SPs in metals; however, it should be understood that, in 
contrast to SPs, here we deal with 2D collective excitations, i.e. plasmons. We 
note that even though field profiles of plasmons in graphene and SPs in metals 
look the same, these two systems are qualitatively different since electrons in 
graphene are essentially frozen in the transverse dimension [3H] . This fact and 
the differences in electronic dispersions (linear Dirac cones vs. usual parabolic) 
lead to qualitatively different dispersions of TM modes in these two systems 
[see Fig. 13.1( b) and Fig. 13.2( d)], To find dispersion of plasmons in graphene 
we need the conductivity of graphene a(u,q), which we now proceed to ana- 
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lyze by employing the semiclassical model (38] (in subsection 13. 2. 1 j) . RPA and 
number conserving relaxation-time approximation [31] (in subsection I3.2.2p . 
and by estimating the relaxation-time due to the influence of electron-phonon 
coupling [35] on the optical conductivity [36] (in subsection 13. 2.3[) . 



3.2.1 Semiclassical model 

For the sake of the clarity of the presentation, we first note that by employing 
a simple semi-classical model for the conductivity (see Ref. [38]), one obtains 
a Drude-like expression: 

< \ i (Q K\ 

■Knr u + it 1 

(the semiclassical conductivity does not depend on q). Here r denotes the 
relaxation-time (RT), which in a phenomenological way takes into account 
losses due to electron-impurity, electron-defect, and electron-phonon scatter- 
ing. Equation (13. 5p is obtained by assuming zero temperature T « 0, which 
is a good approximation for highly doped graphene considered here, since 
Ep ^> ksT. From Eqs. (13. 4p and ( 13. 5 p it is straightforward to obtain plasmon 
dispersion relation: 

q(u}) = — (1 + — )u , (3.6) 

e z Lp too 

as well as losses, 

Mq Itxct 
^sq A a 



Viir 



In order to quantify losses one should estimate the relaxation time r. If the 
frequency u is below the interband threshold frequency uJi nter , and if u < ouo P h, 
then both interband damping and plasmon decay via excitation of optical 
phonons together with an electron-hole pair are inactive. In this case, the 
relaxation time can be estimated from DC measurements [5], i.e., it can be 
identified with DC relaxation time which arises mainly from impurities (see 
Refs. [5]). It is reasonable to expect that impurity related relaxation time will 
not display large frequency dependence. In order to gain insight into the losses 
by using this line of reasoning let us assume that the doping level is given by 
Ep = 0.64 eV (corresponding to electron concentration of n = 3 x 10 13 cm -2 ); 
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the relaxation time corresponds to DC mobility fi = 10000 cm 2 /Vs measured in 
Ref. [5]: Tbc = iih^frvn j ev p = 6.4 x 10 _13 s. As an example, for the frequency 
hu = 0.155 eV (Xair = 8/im), the semiclassical model yields 9ftg/3g ~ 151 for 
losses and X a i r / A p ~ 42 for wave localization. Note that both of these numbers 
are quite favorable compared to conventional SPs [e.g., see Fig. 13.1( c)], It will 
be shown in the sequel that for the doping value Ep = 0.64 eV this frequency 
is below the interband loss threshold, and it is evidently also smaller than the 
optical phonon loss threshold fiujoph ~ 0.2 eV, so both of these loss mechanisms 
can indeed be neglected. 



3.2.2 RPA and relaxation-time approximation 

In order to take the interband losses into account, we use the self-consistent 
linear response theory, also known as the random-phase approximation (RPA) 
[38] . together with the relaxation-time (finite r) approximation introduced by 
Mermin [31]. Both of these approaches, that is, the collisionless RPA (r — > oo) 
[4"U| |4"T) . and the RPA-RT approximation (finite r) [16], have been applied to 
study graphene. In the r — )• oo case, the RPA 2D polarizability of graphene is 
given by [4Tj : 

X(g,^) = ^n(g, w ), (3.8) 



where 



SI J— 1 ?W + J5ni,k- #na,k+q + M7 
k,ni,ri2 

x |(m,k| e -^ r |n 2 ,k + q)| 2 . (3.9) 

Here f(E) = ( e ( E ~ E F)/kBT _|_ j g ^ e p erm j distribution function, Ep is the 
Fermi energy and factor 4 stands for 2 spin and 2 valley degeneracies. Note 
that polarizability x(g, u) is simply related to the density- density response 
function x(g,o;), introduced in chapter El since Ii{q,u) = —x(q,oo). 

Now, in Eq. (13.81) u is given an infinitesimally small imaginary part which 
leads to the famous Landau damping; that is, plasmons can decay by exciting 
an electron-hole pair (interband and intraband scattering) as illustrated in 
Fig. 13.2( c). The effects of other types of scattering (impurities, phonons) can 
be accounted for by using the relaxation-time r as a parameter within the 
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RPA-RT approach [M], which takes into account conservation of local electron 
number. Within this approximation the 2D polarizability is 



-/ \ 1 + i/ur )x {q,u + it) 

XAq,u) = 1 , /■/ — w , ~, ... v 3.10 



The 2D dielectric function and conductivity are respectively given by (see 
and 

a RPA (q,uj) = -iux T (q,ui). (3.11) 

We note here that throughout the text only 7r-bands are taken into consid- 
eration; it is known that in graphite, higher a-bands give rise to a small 
background dielectric constant [13] at low energies, which is straightforward 
to implement in the formalism. Using Eqs. (13. 4 p and (13. lip we obtain that 
the properties of plasmons (i.e., dispersion, wave localization and losses) can 
be calculated by solving 

e RPA (q,u) = 0, (3.12) 

with complex wave vector q = qi + iq 2 . The calculation is simplified by lin- 
earizing Eq. (I3.12p in terms of small q2/qi, to obtain, 



2 



■M[n( qi ,u)} = 0, (3.13) 



2 2e gi 
for the plasmon dispersion, and 

gjn^oQ] + T £^[U(q u u)] + ±&[n( qi ,u)(i - n( gi ,o;))/n( gi ,o)] 

q2 ~ j-n[n(qi^)} - im(qi^)} 

(3.14) 

yielding losses. Note that in the lowest order the dispersion relation (and 
consequently \ a ir/\ and the group velocity v g ) does not depend on r. This 
linearization is valid when q 2 q\] as the plasmon losses increase, e.g., after 
entering the interband regime [the rose area in Fig. I3.2f d)]. results from Eqs. 
(I3.13P and (I3.14p should be regarded as only qualitative. The characteristic 
shape of the plasmon dispersion is shown in Fig. 13.2( d). Note that the semi- 
classical model and the RPA model agree well if the system is sufficiently 
below the interband threshold [for small q, u)(q) ~ y/q as in Eq. (13. 6p ], By 
comparing Figs. 13.2( d) and 13.1( b) we see that the dispersion for SPs on silver- 
dielectric surface qualitatively differs from the plasmon dispersion in graphene 
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[39] . While SPs' dispersion relation approaches an asymptote (u —> usp) 
for large q values [Eq. (13. ip ]. graphene plasmon relation gives u)(q) which 
continuously increases [Fig. 13.2( d)]. 

Theoretically predicted plasmon losses $lq/^sq and wave localization X a i r /\ 
are illustrated in Fig. 13.31 for doping level Ep = 0.135 eV and relaxation 
time r = 1.35 x 1CT 13 s. We observe that for this particular doping level, for 
wavelengths smaller than Xinter ~ 7.7 /xm, the system is in the regime of high 
interband losses (rose shaded region). Below the interband threshold, both 
losses and wave localization obtained by employing RPA-RT approach are 
quite well described by the previously obtained semiclassical formulae. Since 
the frequencies below the interband threshold are (for the assumed doping 
level) also below the optical phonon frequency, the relaxation time can be 
estimated from DC measurements. 

(a) frequency [THz] (b) frequency [THz] ( c ) frequency [THz] 
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Figure 3.3: Properties of plasmons in doped graphene. Solid-lines are obtained 
with the number- conserving RPA calculation, and the dashed lines with the 
semiclassical approach. Losses (a), field localization (wave "shrinkage") (b), 
and group velocity (c) for doping E F = 0.135 eV, and relaxation time r = 
1.35 x 10~ 13 s, which corresponds to the mobility of 10000 cm 2 /Vs. The upper 
scale in all figures is frequency v = u/2n, whereas the rose shaded areas denote 
the region of high interband losses. 

At this point we also note that in all our calculations we have neglected 
the finite temperature effects, i.e., T « 0. To justify this, we note that for 
doping values utilized in this paper the Fermi energies are 0.135 eV~ 5.2A;sT r 
(n = 1.35 x 10 12 cm" 2 ) and 0.64 eV^ 2hk B T r (n = 3 x 10 13 cm" 2 ) for room 
temperature T r = 300 K. The effect of finite temperature is to slightly smear 
the sharpness of the interband threshold, but only in the vicinity (~ ksT r ) of 
the threshold. 

By increasing the doping, Ep increases, and the region of interband plas- 
monic losses moves towards higher frequencies (smaller wavelengths). How- 
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ever, by increasing the doping, the interband threshold frequency will eventu- 
ally become larger than graphene's optical phonon frequency lOoph- there will 
exist an interval of frequencies, uo P h < us < Winter, where it is kinematically 
possible for the photon of frequency u to excite an electron-hole pair together 
with emission of an optical phonon. This second order process can reduce the 
relaxation time estimated from DC measurements and should be taken into 
account, as we show in the following subsection. 



3.2.3 Losses due to optical phonons 

In what follows, we estimate and discuss the relaxation time due to the electron- 
phonon coupling. This can be done by using the Kubo formula which has 
been utilized in Ref. [36] to calculate the real part of the optical conductivity, 
$lcr(uj,q = 0). The calculation of conductivity 9?cr(u;,0) involves the electron 
self-energy whose imaginary part expresses the width of a state with 

energy E, whereas the real part corresponds to the energy shift. Let us as- 
sume that the electron self-energy stems from the electron-phonon coupling 
and impurities, 

E{E) = E e „ ph {E) + E imp (E). (3.15) 

For Ee-ph we utilize a simple yet fairly accurate model derived in Ref. [35]: If 

\E — Ep\ > frujoph, then 

3£ e _ p/l (£) = i\E - sgn(£ - E F )hw 0p h\, (3.16) 

while elsewhere QH e - p h{E) = 0; the dimensionless constant 7 = 18.3 x 10 -3 
[35] is proportional to the square of the electron-phonon matrix element [35] , 
i.e., the electron-phonon coupling coefficient. In order to mimic impurities, 
we will assume that Q : Si mp (i?) is a constant (whose value can be estimated 
from DC measurements). The real parts of the self-energies are calculated 
by employing the Kramers- Kronig relations. In all our calculations the cut- 
off energy is taken to be 8.4 eV, which corresponds to the cut-off wavevector 
k c = it /a, where a = 2.46 A. By employing these self-energies we calculate 
the conductivity §t,a(u, q = 0), from which we estimate the relaxation time by 
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using Eq. (13.51) . i.e., 



r(w) 



e 2 Ei 



(3.17) 



7rft 2 u; 2 Kcx(u;, 0) 

for the region below the interband threshold; in deriving (13.171) we have as- 
sumed ra>l. 

Figure 13.41 plots the real part of the conductivity and the relaxation time 
for two values of doping: Ep = 0.135 eV {n = 1.35 x 10 12 cm -2 , solid line) 
and Ep = 0.64 eV (n = 3 x 10 13 cm -2 , dashed line). In order to isolate 
the influence of the electron-phonon coupling on the conductivity and plas- 
mon losses, the contribution from impurities is assumed to be very small: 
'^s'Ei m p(E) = 10~ 6 eV. The real part of the conductivity has a universal value 
(To — 7re 2 /2/i above the interband threshold value hw = 2Ep (for q = 0), 
e.g., see [37J HI]. We clearly see that the relaxation time is not affected by 
the electron-phonon coupling for frequencies below uo P h, that is, we conclude 
that scattering from impurities and defects is a dominant decay mechanism for 
uj < ujoph (assuming we operate below the interband threshold). However, for 
u > ujoph-, the relaxation times in Fig. 13.41 are on the order of 10~ 14 — 10~ 13 s, 
indicating that optical phonons are an important decay mechanism. 



1 1 

b 0.5 
Pi 





(a) 


















Figure 3.4: (a) The real part of the conductivity in units of ctq = ire 2 /2h in 
dependence of frequency hu/Ep, and (b) the corresponding relaxation time as 
a function of wavelength. The contribution to 9fter(o;) from impurities is chosen 
to be negligible. The displayed graphs correspond to two different values of 
doping which yield E F = 0.135 eV (solid blue line), and Ep = 0.640 eV 
(dashed red line). The position of the optical phonon frequency ftuoph ~ 0.2 
eV is depicted by the dotted vertical line in (b); dot-dashed lines depict the 
values of wavelengths corresponding to 2Ep, that is, the interband threshold 
value (for q = 0) for the two doping concentrations. 

It should be emphasized that the exact calculated values should be taken 
with some reservation for the following reason: strictly speaking, one should 
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calculate the relaxation times r(oj,q) along the plasmon dispersion curve given 
by Eq. (13.1 3ft ; namely the matrix elements which enter the calculation depend 
on q, whereas the phase space available for the excitations also differ for q = 
and q > 0. Moreover, the exact value of the matrix element for electron 
phonon coupling is still a matter of debate in the community. Therefore, the 
actual values for plasmon losses could be somewhat different for oj > ojo v h- 
Nevertheless, fairly small values of relaxation times presented in Fig. 13.41 for 
oj > oooph indicate that emission of an optical phonon together with an electron- 
hole pair is an important decay mechanism in this regime. Precise calculations 
for q > and oj > ojo P h are a topic for a future paper. 



(a) frequency [THz] (b) frequency [THz] (c) frequency [THz] 
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Figure 3.5: Properties of plasmons in doped graphene. Solid-lines are obtained 
with the number- conserving RPA calculation, and the dashed lines with the 
semiclassical approach. Losses (a), field localization (wave "shrinkage") (b), 
and group velocity (c) for doping E F = 0.64 eV; losses are calculated by 
using the relaxation time r -1 = r^ c + T~_} ph , where r DC = 6.4 x 1CT 13 s, 
and T e -ph is the relaxation time from the electron-phonon coupling for the 
given parameters. In the white regions (right regions in all panels), losses 
are determined by Tdc- in the yellow shaded regions (central regions in all 
panels), losses are determined by the optical phonon emission, i.e., T e _ p h- The 
rose shaded areas (left region in all panels) denote the region of high interband 
losses. Dotted vertical lines correspond to the optical phonon frequency 0Jo P h ~ 
0.2 eV. The upper scale in all figures is frequency v = oj/2tt. See text for 
details. 

Plasmonic losses and wave localization calculated from the RPA-RT approx- 
imation are illustrated in Fig. 13.51 for doping level Ep = 0.64 eV and the 
relaxation time r given by r _1 = + t~_ h , where r DC = 6.4 x 10 -13 s 
(mobility 10000 cm 2 /Vs), whereas T e _ p h is frequency dependent and corre- 
sponds to electron-phonon coupling assuming very clean samples [see dashed 
line in Fig. 13.4( b)]. Interband losses [left (rose shaded) regions in all panels] 
are active for wavelengths smaller than Aj nter ~ 1.7 fiva. In the frequency in- 
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terval Ui nte r > uj > uJo P h [central (yellow shaded) regions in all panels], the 
decay mechanism via electron phonon coupling determines the loss rate, i.e., 
t m r e ^ p h- For uj < uo P h [right (white) regions in all panels], the DC relaxation 
time tdc can be used to estimate plasmon losses. 

It should be noted that the mobility of 10000 cm 2 /Vs could be improved, 
likely even up to mobility 100000 cm 2 /Vs [H], thereby further improving plas- 
mon propagation lengths for frequencies below the optical phonon frequency. 
However, for these larger mobilities the calculation of losses should also include 
in more details the frequency dependent contribution to the relaxation time 
from acoustic phonons (this decay channel is open at all frequencies); such a 
calculation would not affect losses for ui > uo P h where optical phonons are 
dominant. 

3.3 Conclusion and Outlook 

In conclusion, we have used RPA and number-conserving relaxation-time ap- 
proximation with experimentally available input parameters, and theoretical 
estimates for the relaxation-time utilizing electron-phonon coupling, to study 
plasmons and their losses in doped graphene. We have shown that for suf- 
ficiently large doping values high wave localization and low losses are simul- 
taneously possible for frequencies below that of the optical phonon branch 
uj < ujoph (i-e., Epi asmon < 0.2 eV). For sufficiently large doping values, there is 
an interval of frequencies above Uoph and below interband threshold, where an 
important decay mechanism for plasmons is excitation of an electron-hole pair 
together with an optical phonon (for uj < u>o P h this decay channel is inactive); 
the relaxation times for this channel were estimated and discussed. We point 
out that further more precise calculations of plasmon relaxation times should 
include coupling to the substrate (e.g., coupling to surface-plasmon polaritons 
of the substrate), a more precise shape of the phonon dispersion curves, and 
dependence of the relaxation time via electron-phonon coupling on q > (see 
subsection I3.2.3P . 

The main results, shown in Figures 13.31 and 13.51 point out some intriguing 
opportunities offered by plasmons in graphene for the field of nano-photonics 
and metamaterials in infrared (i.e. for uj < uioph)- For example, we can see in 
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those figures that high field localization and enhancement \ a ir/\ ~ 200 [see 
Figure IBTBT b)] are possible (resulting in \ < 50 nm), while plasmons of this 
kind could have propagation loss-lengths as long as ~ 10A P [see Fig. 13.5( a)]: 
these values (albeit at different frequencies) are substantially more favorable 
than the corresponding values for conventional SPs, for example, for SPs at 
the Ag/Si interface \ a i r /\ ~ 20, whereas propagation lengths are only ~ 
0.1A sp [see Fig. 13.1( c)], Another interesting feature of plasmons in graphene is 
that, similar to usual SP-systems [TSj, wave localization is followed by a group 
velocity decrease; the group velocities can be of the order v g = 10~ 3 — 10~ 2 c, 
and the group velocity can be low over a wide frequency range, as depicted 
in Figs. 13.3( c) and I3.5f c). This is of interest for possible implementation of 
novel nonlinear optical devices in graphene, since it is known that small group 
velocities can lead to savings in both the device length and the operational 
power [15]; the latter would also be reduced because of the large transversal 
field localization of the plasmon modes. 
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Chapter 4 



Transverse-electric plasmons 



In Chapter [3] we were studying longitudinal charge density oscillations i.e. 
longitudinal plasmons or TM modes. However, due to unusual electron dis- 
persion, graphene can also support transverse plasmons or TE modes |^8j . 
These excitation are possible only if the imaginary part of the conductivity of 
a thin sheet of material is negative [IB]- On the other hand, such a conduc- 
tivity requires some complexity of the band structure of the material involved. 
For example, TE plasmons cannot occur if the 2D material possesses a single 
parabolic electron band. From this perspective, bilayer graphene, with its rich 
band structure and optical conductivity (e.g., see [19] and references therein), 
seems as a promising material for exploring the possibility of existence of TE 
plasmons. Here we predict the existence of TE plasmons in bilayer graphene. 
We find that their plasmonic properties are much more pronounced in bilayer 
than in monolayer graphene, in a sense that the wavelength of TE plasmons 
in bilayer can be smaller than in monolayer graphene at the same frequency. 

Throughout this work we consider bilayer graphene as an infinitely thin sheet 
of material with conductivity <r(q, u). We assume that air with e r = 1 is above 
and below bilayer graphene. Given the conductivity, by employing classical 
electrodynamics, one finds that self-sustained oscillations of the charge occur 
when (see [IB] and references therein) 
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for TM modes, and 

1 — &0^1 =n (4.2) 

2^ - u 2 /c 2 

for TE modes. The TM plasmons can considerably depart from the light line, 
that is, their wavelength can be considerably smaller than that of light at the 
same frequency. For this reason, when calculating TM plasmons it is desirable 
to know the conductivity as a function of both frequency u and wavevector q. 
However, it turns out that the TE plasmons (both in monolayer [H] and bilayer 
graphene, as will be shown below) are quite close to the light line q = co/c, 
and therefore it is a good approximation to use <j(u) = a(q = 0, u). Moreover, 
these plasmons are expected to show strong polariton character, i.e., creation 
of hybrid plasmon-photon excitations. At this point it is worthy to note that if 
the relative permittivity of dielectrics above and below graphene are sufficiently 
different, so that light lines differ substantially, then TE plasmon will not exist 
(perhaps they could exist as leaky modes). 



4.1 Optical conductivity of bilayer graphene 

The conductivity cr(uj) = 9?cr(u;) + i^sa(u) is complex, and plasmon disper- 
sion is characterized by the imaginary part Qa(uj), whereas 9fax(a;) determines 
plasmon losses, or more generally absorption of the sheet. From Eq. (14 .2p it 
follows that the TE plasmons exist only if Qa(u) < [4"8] . 




Figure 4.1: The band-structure of bilayer graphene. The two upper bands (as 
well as the two lower bands) are perfectly nested and separated by 7 ~ 0.4 eV; 
<?o = j/hvF- Horizontal line depicts one possible value of the Fermi level, and 
arrows denote some of the possible interband electronic transitions. See text 
for details. 
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In order to calculate the imaginary part of the conductivity, we employ 
Kramers-Kronig relations and the calculation of absorption by Nicol and Car- 
botte [19], where Ua(ou) [see Eqs. (19)-(21) in Ref. [19]] was calculated by 
using the Kubo formula. The optical conductivity has rich structure due to 
the fact that the single-particle spectrum of graphene is organized in four bands 



where vp = 10 6 m/s, the parameter 7 « 0.4 eV is equal to the separation 
between the two conduction bands (which is equal to the separation between 
the valence bands). The band structure (14. 3 p is calculated from the tight 
binding approach, where Vp is connected to the nearest-neighbour hopping 
terms for electrons to move in each of the two graphene planes, and the distance 
between Carbon atoms in one monolayer (see Ref. [49J), whereas 7 is the 
hopping parameter corresponding to electrons hoping from one layer to the 
other and vice versa [19]. The two graphene layers are stacked one above the 
other according to the so-called Bernal-type stacking (e.g., see Ref. [50]). We 
emphasize that the perfect nesting of bands gives rise to the stronger plasmon 
like features of TE plasmons in bilayer than in monolayer graphene. The four 
bands are illustrated in Fig. 14.11 along with some of the electronic transitions 
which result in absorption. Absorption depends on 7 and the Fermi level \x\ 
the latter can be changed by applying external bias voltage. 



The imaginary part of the conductivity can be calculated from dta(oj) by 
using the Kramers-Kronig relations 



given by [4"9| . 




(4.3) 




(4.4) 
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which yields 
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g(x,y,z) 
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(4.6) 



o"o = e 2 /2h, Q(x) — 1 if x > and zero otherwise, and f2 = faj. Here we 
assume zero temperature T f» 0, which is a good approximation for sufficiently 
doped bilayer graphene where \i 3> /cbT. Formulae (14. 5p and (14. 6 p are used to 
describe the properties of TE plasmons. 

In Figure 14.21 we show the real and imaginary part of the conductivity for 
two different values of the Fermi level: /t = O.47 and /t = O.97 (we focus on 
the electron doped system \x > 0). Because plasmons are strongly damped by 
interband transitions, it is instructive at this point to discuss the kinematical 
requirements for the excitation of electron-hole pairs. If the doping is such 
that \x < 7/2, a quantum of energy ftw (plasmon or photon) with in-plane 
momentum q = can excite an electron- hole pair only if fkv > 2/i (excitations 
from the upper valence to the lower conduction band shown as red dot-dashed 
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line in Fig. 14. ip . If \i > 7/2, the (q = 0, w)-quantum can excite an electron-hole 
pair only for huj > 7 (excitations from the lower to the upper conduction band 
shown as green solid lines in Fig. 14. II occur at fiuj = 7). If the plasmon/photon 
has in-plane momentum q larger than zero, then interband transitions are 
possible for smaller frequencies (see blue dashed lines in Fig. I4.ip . There is a 
region in the (q, o;)-plane where electron-hole excitations are forbidden due to 
the Pauli principle. Because plasmons are strongly damped by these interband 
transitions (this is Landau damping), in our search for the TE plasmons, we 
focus on their dispersion curve in the regime where electron-hole pair formation 
is inadmissible (via first-order transition). 



4.2 Transverse-electric plasmon dispersion in 
bi layer graphene 




Figure 4.2: The real (red dotted lines) and imaginary (blue solid lines) part 
of the conductivity of bilayer graphene for two values of doping: \i = O.47 
(a), and /i = O.97 (b). The conductivity is in units of er = e 2 /2h, and 
the frequency is in units of Uq — j/H. The 5-functions in 3ta(u) at uj = 
(intraband transitions) and oj — j/H (transitions from the lower to the upper 
conduction band depicted as green solid arrows in Fig. 14. 1 j) are not shown (see 

TO- 



In Figure 14.31 we show the plasmon dispersion curves for fi = O.47 and 
[i = O.97; in the spirit of Ref. [3S], we show Aq = q — u/casa function of 
frequency uj. Plasmons are very close to the light line and thus one can to a 
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very good approximation write the dispersion curve as 



Aq 



CO 



(4.7) 



To the left (right) of the vertical red dotted line in Fig. I4.3[ plasmon damping 
via excitation of electron-hole pairs is (is not) forbidden. 
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Figure 4.3: The plasmon dispersion curve Aq = q — u/c vs. oj for \x = 
O.47 (a), and fi = O.97 (b) is shown as blue solid line. To the right of the 
vertical red dotted lines plasmons can be damped via excitation of electron- 
hole pairs, whereas to the left of this line these excitations are forbidden due 
to the Pauli principle. Black dashed line in (b) (which closely follows the blue 
line) corresponds to Eq. (I4.8p . The wave vector is in units of q = •y/hvp, and 
the frequency is in units of ujq = j/h. 

For /1 = O.47, 5s a (u) is smaller than zero for u in an interval of frequencies 
just below 2/i. From the leading term in 5sa(u) we find that departure of 
the dispersion curve from the light line is logarithmically slow: Ago<^< 7 /2 oc 
[log \huj — 2/x|] 2 . The same type of behavior occurs in monolayer graphene |4"8"] . 

However, for fi = O.97, one can see the advantage of bilayer over monolayer 
graphene in the context of TE plasmons. The conductivity 5sa(uj) is smaller 
than zero in an interval of frequencies below 7. In this interval, the most 
dominant term to the conductivity is the last one from Eq. f)4.5p . that is, 
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2^2 „3 



r 



(4.8) 



This approximation is illustrated with black dashed line in Fig. I4.3[ and it 
almost perfectly matches the dispersion curve. Note that the singularity in 
5sa(uj) at hbj = 7 is of the form 1/(7 — hw), whereas the singularity at fojj — 1\i 
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is logarithmic (as in monolayer graphene [H]). As a consequence, the departure 
of the dispersion curve from the light line in bilayer graphene is much faster for 
fi > 7/2 than for fi < 7/2, and it is faster than in monolayer graphene as well 
[note the two orders of magnitude difference between the abscissa scales in Figs. 
14.3( a) and (b)]. Thus, we conclude that more pronounced plasmonic features 
of TE plasmons (shrinking of wave length which is measured as departure of 
q from the light line) can be obtained in bilayer graphene. The term in $sa(uj) 
which is responsible for TE plasmons for \i > 7/2 corresponds (via Kramers- 
Kronig relations) to the absorption term b(fi)5(hu — 7) [H], which arises from 
the transitions from the first to the second valence band (shown as green solid 
arrows in Fig. 14. ip . which are perfectly nested and separated by 7. Thus, 
this unique feature of bilayer graphene gives rise to TE plasmons with more 
pronounced plasmon like features than in monolayer graphene. 

Before closing this chapter, let us discuss some properties and possible ob- 
servation of TE plasmons. First, note that since the electric field oscillations 
are both perpendicular to the propagation vector q, and lie in the bilayer 
graphene plane, the electric current j = cr(u;)E is also perpendicular to q. 
Thus, j • q = 0, and the equation of continuity yields that the charge density is 
zero (i.e., one has self-sustained oscillations of the current). In order to excite 
plasmons of frequency u with light of the same frequency, one has to somehow 
account for the conservation of the momentum which is larger for plasmons. 
Since the momentum mismatch is relatively small, the standard plasmon ex- 
citation schemes such as the prism or grating coupling methods (e.g., see [12] 
and references therein) could be used for the excitation of these plasmons. 

To conclude this chapter, we have predicted the existence of transverse elec- 
tric (TE) plasmons in bilayer graphene. Since they exist very close to the light 
line, these plasmons are expected to show strong polariton character, i.e., mix- 
ing with photon modes. However, due to the perfectly nested valence bands of 
bilayer graphene, their dispersion departs much more from the light line than 
in monolayer graphene. 



59 



60 



Chapter 5 

Plasmon-phonon coupling 



In this chapter we analyze the coupling of plasmons with intrinsic optical 
phonons in graphene by using the self-consistent linear response formalism. 
We find that longitudinal plasmons (LP) couple only to transverse optical 
(TO) phonons, while transverse plasmons (TP) couple only to longitudinal 
optical (LO) phonons. The LP-TO coupling is stronger for larger concentration 
of carriers, in contrast to the TP-LO coupling (which is fairly weak). The 
former could be measured via current experimental techniques. Thus, plasmon- 
phonon resonance could serve as a magnifier for exploring the electron-phonon 
interaction, and for novel electronic control (by externally applied voltage) 
over crystal lattice vibrations in graphene. 

To analyze plasmon-phonon coupling let us start with the Hamiltonian for 
the Dirac electrons in graphene 

H e = tivpcr ■ k, (5.1) 

where vp = 10 6 m/s, k = (k x ,k y ) = — iV is the wave-vector operator, a = 
(a x ,a y ), and a x>y are the Pauli spin matrices. We label the eigenstates of 
Hamiltonian H e by |s,k) and the appropriate eigenvalues by E s ^, = sHvf\^\, 
where s = 1 for the conduction band and s = — 1 for the valence band. 

The long-wavelength in-plane optical phonon branch in graphene consists of 
two modes (LO and TO) which are effectively dispersionless and degenerate at 
energy hwo = 0.196e^. Let u(R) = [u^R) — ub(R)]/\/2 denote the relative 
displacements of the sub-lattice atoms A and B of a unit cell specified by a 
coordinate R [see Fig. 15.1( c)], Then, in the long- wavelength limit R can be 
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replaced by a continuous coordinate r and we have 



u « = E 7^ g ^ eiqr ' (5 - 2) 



where N is the number of unit cells, M is the carbon atom mass, q = 
g(cos0 q , sin0 q ) is the phonon wave vector, \x = L,T stands for the polar- 
ization, and the polarization unit vectors are e^ q = i(cos0 q , sin0 q ), and 
e Tq = «(— sin q , cos q ). The displacement vector u(r) is parallel (perpen- 
dicular) to the phonon propagation wave vector q for LO (TO, respectively) 
phonons [see Fig. 15.1( c)]. The phonon Hamiltonian is given by 

R ph = \ E( P ^ P ^ + ^oQUQ m ), (5-3) 

where Q m and P m denote phonon coordinate and momentum. The electron- 
phonon interaction takes a peculiar form in graphene (see chapter [2]): 

H e . ph = -V2^<t x u(r), (5.4) 

where a x u = a x u y — a y u x , b = 0.142 nm is the nearest carbon atoms distance, 
and (3 = 2. We find it convenient to write Eq. (15 .4p as 

H e „ ph = L 2 F^j q x e m Q m (5.5) 
m 

where j q = —evpL~' 1 ae^ % ° Lr is the single-particle current- density operator, L 2 
is the area of the system, e is charge of the electron, and F = -^=. 

J ' to ' eb 2 VNM 

The electromagnetic field in the plane of graphene is completely described by 
the vector potential A = J2 m e Mq^Mq e * qr (scalar potential is gauged to zero, 
time dependence is implicitly assumed, and /i = L,T denote polarizations). 
The interaction with Dirac electrons is obtained by substitution Kk — > Kk + eA 
in Eq. ( 15. ip . which leads to 

He-ern = eVpCT ■ A = -1? ^ Jq ' ^q^W ( 5 " 6 ) 

By comparing Eqs. (15 .4p and (I5.6P it follows that electron-phonon interaction 
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Figure 5.1: (a) Schematic illustration of the lattice structure with two sublat- 
tices (A and B). (b) The two degenerate Dirac cones are centered at K and K' 
points at the edge of the Brillouin zone, (c) A displacement of lattice atoms 
u(r) is parallel (perpendicular) to the propagation wave vector q of a LO (TO) 
phonon. (d) The displacement u(r) creates an effective vector potential A e ff 
perpendicular to u(r) (the sign of A e g for the K' point is opposite to that for 
the K point). 

can be regarded as a presence of an effective vector potential 

A cff = F ^(e Tq g Lq - e Lq Q Tq )e iqr , (5.7) 
q 

that is, H e _ p h = evpcr- A e g. It is evident that A c g-u(r) = that is the effective 
vector potential A e g is perpendicular to u(r) as illustrated in Figs. 15. lif e) and 
(d) (see also Ref. [50]), which is responsible for the mixing of polarizations in 
plasmon-phonon coupling. 

As a first pass, let us ignore the phonons and focus on the Hamiltonian 
H = H e + H e _ em . Without an external perturbation, the electrons in graphene 
fill the Fermi sea according to the Fermi distribution function / s k. A field 
A m {oj) oscillating at frequency to will induce an average current density (up 
to a linear order in the vector potential) 

(J M (q,w)> = -% M (q,o;)A Mq (a/), (5.8) 

where the current-current response function (including 2-spin and 2- valley de- 
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^eneracy) is given by [56 



X/1 (q, co) = 4L 2 z ^/" k 4-- 

x |( Sl k|j q -e; q | S2 k + q)| 2 . (5.9) 

For the response function x M (q, o>) we utilize the analytical expression from 
Ref. [54j. The subtlety involved with the divergence in Eq. (15. 9p is solved by 
subtracting from Xl(q, [xr(q, ^)] the value Xz^q? w = 0) [xr(q - > 0, cj — 0)] 
to take into account that there is no current response to the longitudinal 
[transverse] time [time and space] independent vector potential, see [5H [55] 
for details. We would like to note that when working with the current-current 
response function, rather than with the density-density response function, the 
nature of the plasmon-phonon interaction (especially the mixing of polariza- 
tions as shown below) is far more transparent. 

Next, it is straightforward to show from the Maxwell equations that an 
electric current oscillating in a two-dimensional plane will induce a vector 
potential 

(An(u>)) = (^(q^))^^^, (5.10) 



and 



(A Tq H) = (J r (q, (5.11) 

2y q 2 — u 2 /c 2 



where we have assumed that graphene is suspended in air and that there are 
no other sources present in space. This induced vector potential in turn acts 
on electrons in graphene through the interaction Hamiltonian H e _ em which 
can result in plasmons - self-sustained collective oscillations of electrons. From 
Eqs. (15. 8p and (15. 10p we get the dispersion relation for longitudinal plasmons 

HUE] 



1 " o..2 X*(q, w) = 0. (5.12) 



\/ q 2 — i0 2 /c 2 
2cu 2 e 

From Eqs. ( 15. Sp and (15. lip we get the dispersion relation for transverse plas- 
mons [H] 

Xr(q,") = 0. (5.13) 



2vV - u 2 /c 2 

Longitudinal plasmons are also referred to as transverse magnetic modes since 
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they are accompanied by a longitudinal electric (E) and a transverse magnetic 
field (B) in the plane of graphene. Likewise transverse plasmons or transverse 
electric modes are accompanied by a transverse electric and a longitudinal 
magnetic field [48]. Dispersion relation of LP (TP) modes is shown by the 
blue dashed line in Fig. 15.21 (Fig. 15.31 respectively). Finally we note that 
we are primarily interested in non-radiative modes (q > oj/c) in which case 
fields are localized near the graphene plane (z = 0) and decay exponentially: 
E(z),B(z) oc e-HvV-" 2 /" 2 . 

In order to find the plasmon-phonon coupled excitations we consider the 
complete Hamiltonian H = H e + ff e _ em + H e _ p h + H P h- We assume that the 
hybrid plasmon phonon mode oscillates at some frequency uj with wavevec- 
tor q (which are to be found). From the equation of motion for the phonon 
amplitudes Q m one finds [56] 

(io 2 -u 2 )(Q Tn ) = L 2 F(J L (q,u)), (5.14) 

and 

(u 2 -co 2 )(Q L(i ) = -L 2 F(J T (q,co)). (5.15) 

The electron phonon interaction (15. 5 j) is included as the effective vector poten- 
tial (15 .7]) in Eq. (15. SJ) , which from Eq. (15. 8p immediately yields 

(J L (q,u)) = XL (ci,ou)(-(A L(l (u)) + F(Q Tq )), (5.16) 

and 

(Jr(q,w))=XT(q,w)(-<A Tq (a;))-F<Q Lq )). (5.17) 

From Eqs. (I5.14p - (15.171) it is clear that transverse (longitudinal) phonons 
couple only to longitudinal (transverse) plasmons. Apparently, this follows 
from the fact that LO (TO, respectively) phonons are equivalent to oscillations 
of an effective vector potential A c fj [see Eq. (j5.7|) ]. and therefore an effective 
electric field, perpendicular (parallel, respectively) to q. 

Finally using Eqs. (15.10]) . (I5.14p . and ( I5.16P we get the dispersion relation 
for the LP-TO coupled mode 

2 2 L 2 F 2 X L(q,u) , . 

v -u = , (5.18) 
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Figure 5.2: Dispersion lines of hybrid LP-TO plasmon-phonon modes (solid 
lines) and of the uncoupled modes (dashed lines) for two values of doping: (a) 
n = 5 x 10 12 cm -2 , and (b) n = 5 x 10 13 cm -2 . The hybridization is stronger for 
larger doping values. Grey areas denote the region of single-particle damping. 



and from Eqs. (15. lip . (15.151) . and (15.171) dispersion relation for the TP-LO 
coupled mode 

2 2 L 2 F 2 X T{q,u) , . 

u _ w o = 1 iTn 1 (5.19) 



O I 2 2,2 ^^" 

The plasmon dispersions relations (15 . 12[) and f 1 5 . 1 3 j) appear as poles in the 
Eqs. (I5.18P and (I5.19P for the coupled modes, which means that the coupling 
is greatest at the resonance point where plasmon momentum and energy match 
that of the appropriate phonon mode. We denote this point (where the uncou- 
pled plasmon and phonon dispersion cross) by (q c ,uo). One can quantify the 
strength of the coupling effect by calculating the frequency difference between 
the hybrid modes at the wavevector q c in units of the uncoupled frequency 
value: Au/uq. Finally by doping one can change plasmon dispersion which in 
turn changes q c and the strength of the plasmon-phonon coupling. 



The dispersion lines for the hybrid LP-TO modes are shown in Fig. 15.21 
for two values of doping, (a) n = 5 x 10 12 cm -2 , E F = 0.261 eV, k F = 
3.96 x 10 8 nr 1 , and (b) n = 5 x 10 13 cm" 2 , E F = 0.825 eV, k F = 1.25 x 10 9 m" 1 . 
The strength of the coupling increases with increasing values of doping, and 
one has for the case (a) Acu/ujq = 7.5%, and (b) Auj/u = 15.5%. To describe 
graphene sitting on a substrate (say SiC, which is a polar material), one only 
needs to include the dielectric function of the substrate into our calculation. 
In that case plasmons can also couple to surface phonon modes of the polar 
substrate [57]. However, since these surface phonons have sufficiently smaller 
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Figure 5.3: Dispersion lines of hybrid TP-LO plasmon-phonon modes (solid 
lines) and of the uncoupled modes (dashed lines) for two values of doping: 

The plasmon-like 
the ordinate shows 



n 



(a) n = 7.5 x 10 cm -2 , and (b) 
dispersion is very close to the light line q 
Aq = q — uj/c. 



9.5 x 10 11 cm" 2 
u/c; therefore. 



energies than optical phonons in graphene out results are qualitatively un- 
changed in that case. LP-TO hybrid modes could be measured by observing 
the change in the phonon dispersion with the Neutron Spectroscopy or In- 
elastic X-ray Scattering. Alternatively, one could use grating coupler or Elec- 
tron Energy Loss Spectroscopy to measure the shift in the plasmon energy. 
Our results imply that plasmon-phonon coupling could serve to explore the 
electron-phonon interaction (the frequency shifts at resonance are much larger 
then the G peak shift recently measured by Raman Spectroscopy [32]), and 
that by externally appling voltage one can influence the properties of lattice 
vibrations. 

In spite of the fact that the formal derivation of hybrid TP-LO coupled 
modes is equivalent to the derivation of the LP-TO modes, their properties 
qualitatively differ. First, we note that the dispersion of transverse plasmons 
is extremely close to the light line, and we plot Aq = q — u/c vs. frequency 
bj following Ref. jH]. For this reason, transverse plasmons are expected to 
have strong polariton character and they will be hard to distinguish from free 
photons (also, even a small plasmon linewidth will obscure the distinction). 
Moreover, they do not exist in graphene between two dielectrics with suffi- 
ciently different relative permittivity, where the light lines for the dielectrics 
are separated. Next, transverse plasmons exist only in the frequency inter- 
val 2Ep > fojj > 1. 667 Ep [IS], which means that the LO phonon energy 
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must be in the same interval for the hybridization to occur. Figure 15.31 shows 
the dispersion curves of the hybrid TP-LO modes for two values of doping, 
(a) n = 7.5 x 10 11 cm" 2 , E F = 0.101 eV, k F = 1.53 x 10 8 m" 1 , and (b) 
n = 9.5 x 10 11 cm" 2 , E F = 0.114 eV, k F = 1.73 x 10 8 m _1 . We observe that 
the trend here is opposite to that of the LP-TO coupling, as the strength of the 
coupling decreases with increasing doping; specifically, one has for the case (a) 
Auj/uq = 0.17%, and (b) Au/uq = 0.02%. The maximal coupling occurs when 
2E F is just above Huq, and it is zero when huio = 1.667-Ep. We emphasize that 
the strength of the coupling for TP-LO modes is in general much weaker than 
in LP-TO modes. 

Before closing this chapter, we note another interesting result which is cap- 
tured by our calculations. Equations (I5.18P and (15.191) for shifts in the energies 
of TO and LO modes at q = reduce to 

1 + 2^XL,T(0,W) 



which is identical to the result of Ref. [5T], where the coupling of optical 
phonons to single-particle excitations was studied, appart from the imaginary 
term in the denominator which is zero in [51]. This small but qualitative 
difference is consequence of phonon coupling to the radiative electromagnetic 
modes, which increases the phonon linewidth. For example, for the doping 
values of n = 5 x 10 12 cm -2 , 5 x 10 13 cm -2 , and 5 x 10 14 cm -2 , Eq. f)5.20p 
yields 0.005%, 0.07%, and 0.7%, respectively, for the linewidths, while there is 
no linewidth from single-particle damping at these doping values. This effect 
is qualitatively unchanged for graphene sitting on a substrate and could be 
measured by Raman spectroscopy. Finally, we note an interesting solution 
of Eq. (15.191) (valid for suspended graphene): when the hybrid TP-LO mode 
dispersion crosses the light line it has the same energy as the uncoupled phonon 
mode, i.e., u = ujq. In other words, LO phonon at a wavevector q = u>o/c 
decouples from all (single particle and collective) electron excitations, while 
no such effect exists for the TO phonons. 

In conclusion, we have predicted hybridization of plasmons and intrinsic op- 
tical phonons in graphene using self-consistent linear response theory. To the 
best of our knowledge, this is the first study of such resonance in an isolated 2D 
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material. We found that graphene's unique electron-phonon interaction leads 
to unconventional mixing of plasmon and optical phonon polarizations: lon- 
gitudinal plasmons couple exclusively to transverse optical phonons, whereas 
graphene's transverse plasmons couple to longitudinal optical phonons; this 
contrasts plasmon-phonon coupling in all previously studied systems. The 
strength of the hybridization increases with doping in LP-TO coupled modes, 
while the trend is opposite for TP-LO modes. The LP-TO coupling is much 
stronger than TP-LO coupling, and it could be measured by current experi- 
ments, which would act as a magnifier for exploring the electron-phonon in- 
teraction in graphene. This coupling is an even more striking example of a 
breakdown of Born-Oppenheimer approximation in graphene than the recently 
measured stiffening of the Raman G peak [52] • Moreover, plasmon-phonon in- 
teraction can serve to electronically control the frequencies of lattice vibrations 
in graphene, which could have interesting technological implications. 
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Chapter 6 



Near field heat transfer 

6.1 Near field heat transfer between two graphene 
sheets 



In this chapter we analyze the near field heat transfer between two graphene 
sheets mediated by thermally excited plasmon modes and demonstrate that 
there is a large enhancement of heat transfer compared to the far field black 
body radiation. The system we analyze, shown in figure I6.1[ consists of a 
suspended graphene sheet at temperature T\ emitting to another suspended 
graphene sheet held at room temperature T 2 = 300-fT, and a distance D away. 

In chapter [2] we calculated the expression for the radiative heat exchange 
between two graphene sheets. Total heat transfer H = Hff + H n f can be 
conveniently separated into the contribution from the propagating waves in 
the far field 



H ff = - dwte^TO-e^Ta)] / qdq}^h» ff (q,u;), (6-1) 



Here Q(cu,T) = hujjieP^ — 1) is the Boltzman factor, \i stands for s or p 





and evanescent waves in the near field 




/ tMe^Taj-e^rO] / q^^KM") 



(6.2) 
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polarization and functions and h^f are given by: 
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Figure 6.1: (a) Schematic diagram of the radiation transfer problem: a free 
standing sheet of graphene at temperature 7\ is radiating to another free stand- 
ing graphene sheet at temperature T 2 and distance D away, (b) Schematic dia- 
gram of the field profile for even mode, (c) Odd mode, (d) Contour plot of the 
transfer function h p n f for the case of two graphene sheets at the same chemical 
potentials /ii^ = 0.5eV and same temperatures Ti j2 = 300-fT, separated by a 
distance of D — 10 nm. Dashed line denotes the plasmon dispersion relation 
for a single isolated graphene sheet, while poles of transfer function h p n j show 
dispersion relations of the coupled (even and odd) modes of the two sheets. 

Let us first note that graphene is a poor absorber in the far field, since it 
is only one atom thick. Indeed, it was experimentally demonstrated [37J that 
graphene absorbs only \a\ 2 m 2% of the incident light (see equation (12.771) ). 
Since one can also write |a| 2 = l — \r\ 2 — |t| 2 we can simply neglect the far field 
transfer (see equations (16. ip and (I6.3p ). at least compared to the black body 
case which is characterized by |a#.B| 2 = 100%. On the other hand, we will see 
that near field heat transfer can be significantly greater than the black body 
case, if graphene sheets are sufficiently close to allow the tunneling of surface 
modes (plasmons). 

To analyze the near field heat transfer between two graphene sheets let us 
write the p polarization reflection coefficient (12.1311) for a single sheet as r p = 
(1 — e)/e where e = 1 + 7<r/(2e cj) is the dielectric function of graphene 
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We immediately see that poles of r p are located at the plasmon dispersion 
e = which was derived in chapter [3] but we write it here again for the sake 
of clearance: 

a(uj,T) 



We assumed here 7 = ^u 2 /c 2 — q 2 ~ iq since we have shown that plasmon 
dispersion is mainly located in the non-retarded regime q » oj/c. We have 
already pointed out that strong near field heat transfer requires graphene sheets 
to be very close, which in turn allows coupling of two plasmon modes (see figure 
16. ip . In the case when two graphene sheets have identical parameters (r p = rf) 
this coupling results in two new modes: even mode described by an equation 
r P = e qD ^ anc j an QC ^ moc [ e described by an equation r p = —e qD . Naturally, 
when sheets are sufficiently far apart (D >>) or the wave vector is sufficiently 
large (q >>) so that the coupling becomes irrelevant, these two modes become 
degenerate again, described by a pole in r p . Further on, note that these two 
equations can be joined in a single one: 1 — (r p ) 2 e~ 2qD = 0, which is precisely 
a denominator in equation (16.41) which determines the poles of the function 
h p n j. In other words, these two coupled surface modes, strongly enhance and 
dominate the near field heat transfer. Note however from equation (16. 4p that 
h p n j oc e~ 2qD so that graphene sheets have to be very close to have a significant 
near field heat transfer. In other words plasmon surface modes can act as an 
excellent heat conductors, only the graphene sheets have to be very close to 
allow the coupling of exponentially decaying (E oc e~ qD ) plasmon field. 

Finally note that h p n j oc Sr^Qrf , while Qr p has a pole at the bare plasmon 
dispersion so there will be a competition between this factor and a pole at a 
dispersion of a coupled mode. Therefore the function h p n f will increase with 
increasing wave vector since the coupling between the modes will decrease and 
even/odd mode dispersion (and a corresponding pole) will join that of a bare 
plasmon (and a corresponding pole). However, when these dispersions meet 
(q m 1/D) the function h p n j will start to decrease with the wave vector due 
to exponentially decaying factor h p n ^ oc e~ 2qD . At last note that function h^f 
is multiplied with a Boltzman factor Q(oj,T) which shifts everything to lower 
frequencies so there are several competing effects in action which will be hard 
to disentangle in the end when everything gets integrated over all q, u values. 

The same analysis applies to the s polarization however it is easy to see that 
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it will have a minor contribution to the total near field heat transfer. The 
reason for this is a large difference in the character of the plasmon disper- 
sion relations (compare figure l5\2l and figure l5\3|) . On one hand longitudinal 
plasmons, described by a pole in r p , are located in the non- retarded regime 
(q » co/c) with a large density of states, while transverse plasmons, describe 
by a pole in r s , are located in the strongly retarded regime (q ~ oj/c) with a 
tiny density of states. Since each q value can be thought of as a separate heat 
channel, and if graphene sheets are close enough so that all relevant q modes 
are active, the p polarization will have many more heat channels and dominate 
over the s polarization. 

To model graphene we shall use g-independent conductivity which simplifies 
the mathematical calculations and gives a good order of magnitude on the 
heat transfer (see discussion below). In chapter [2] we showed that the total 
conductivity cr(w) = <t d (lu) + crj(u), can be separated into Drude (intraband) 
and interband part, expressed respectively as (see also [6T]): 



where G(e) = sinh(e//c&T)/(cosh(/i/&; b T) + cosh(e//c fc T)), and /i is the chemical 
potential. Various electron scattering processes are taken into account through 
the relaxation time r. From DC mobility measurements in graphene, one ob- 
tains an order-of-magnitude value of r ~ 10~ 13 s. Now, due to fluctuation- 
dissipation theorem h^f oc Qr p and Qr p oc 9ftcr (see equation (12.1471) ). so we 
have to take particular attention to the origin of dissipation (3fta) in our system. 
At zero temperature the situation is very simple since Drude term (intraband 
contribution) and relaxation time r determines the losses for low frequencies, 
while interband contribution is dominant for frequencies above the interband 
threshold {hw = 2/i). However, at finite temperature, interband processes can 
play a leading role even below the absorption threshold u ~ 2/i, particularly for 
small chemical potential where thermal broadening of interband threshold (on 
the order of few kbT) becomes more significant. While the use of g-independent 
expression for graphene conductivity (16.61) for intraband processes is a good 
approximation, one must take care when applying (16. 6p to interband transi- 




(6.6) 
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tions. Here, the contribution from the finite wave-vector becomes important 
since it broadens the interband threshold from 2/i to 2/i — hqvp. On the other 
hand, this is similar to finite temperature effects which also broaden the in- 
terband threshold, so we do not expect a qualitatively different result with 
q- dependent conductivity. 



Figure 6.2: Contour plot of the near field heat transfer between two graphene 
sheets H™* normalized to the far field heat transfer between two black bodies 
Hg B of the same temperatures, in the log scale. Here T 2 = 300K, /ii >2 = O.leV 
and Ti )2 = l(T 13 ,s. 

To quantify the heat exchange in the near field we plot in figure I6T21 the total 
transfer H n f (16. 2p normalized to the transfer between two black bodies in the 
far field. Since the exponentially decaying Boltzman factor shifts all the con- 
tributions to the lower frequencies we will focus on the small values of chemical 
potential. For \x\^ = O.leV, we observe orders-of-magnitude increase in heat 
exchange particularly at small separations (x 1000 for D = 20nm, 7\ = 800/i), 
but also at separations as large as O.lfim. In general, dependence of transfer 
on separation D is non-uniform and does not seem to yield a simple functional 
dependence on the emitter and absorber temperatures (as is the case for two 
black bodies). This efficient heat exchange between two graphene sheets in the 
near field, together with recently reported advances in hot carrier extraction 
from graphene [62], may offer a potential for a novel, hybrid thermophoto- 
voltaic/thermoelectric solid-state heat-to-electricity conversion device. 
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6.2 Near field thermo-photo-voltaics using graphene 
as a thermal emitter 



Here we show that graphene can be used as a thermal emitter in the near field 
thermo-photo-voltaic (TPV) system resulting in high efficiencies and power 
densities. The near field heat transfer is mediated by thermally excited plas- 
mon modes in graphene similarly to the situation in the last section. 

The system we analyze consists of a hot graphene emitter at a temperature 
T\ and a photo-voltaic (PV) cell held at room temperature T2 = 300K and 
distance D away from graphene. It is interesting to note that the expression 
for the near field heat transfer between graphene and a PV cell is given by the 
same expression (16. 2p for the graphene to graphene heat transfer 



1 f°° 
H nf = — J 



fajJ fad 



POO 

/ qdqh p nf (q,u), (6.7) 
J oj/c 



where 



Q r p Cj r p e - 2 <i D 

h nM> w ) = \i- r Pr\e--* D \* ' (6 ' 8) 
Here we have neglected the contribution from s polarization, and assumed 7 = 
^/w 2 /c 2 — q 2 « iq since the near field heat transfer is mediated by graphene 
plasmon modes in the non- retarded regime (q » uj/c). To model graphene 
we use g-independent conductivity cr(cu) from equation (16.61) as before, and the 
reflection coefficient is given by equation (12 . 1 3 1 [) which we write again for the 
sake of clearance: 

iqa(u)) 

^) = -Mh- (6-9) 

' 2e u 

Above hot graphene emitter we now have PV cell which we model as a simple 
direct band-gap semiconductor with parameters: 

/ a \ 2 , u < u g 

e 2 [u>) = in + i— where a(u) = < rj— (6.10) 

V 2k oJ ao y^f ,u>u g 

Here n is the refractive index, ko = 2ir/\ = c/u is the photon wavelength 
in vacuum, and u g is the bandgap frequency. Specifically we will discuss the 
case of indium antimonide (InSb) with parameters u g = 0.1 7eV and ~ 
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0.7 x 10 4 cm 1 (at room temperature [M])- Finally the reflection coefficient of 
the PV cell, in the non- retarded regime (q » uj/c), is simply given by [HZ] 

r ^ = ^RTT- (6 ' n) 

When the PV cell is biased at a voltage V , we can express the total radiative 
power exchange as [65] 

1 f°° T hw hw 1 f°° 

On the other hand we can also write the total photon flux into the PV cell as 



Jph r, 

7T Z 



poo roo 

/ duJ eftfr, _ i - e W^Vo) _ l I Qdqh p nf (q,u). (6.13) 

JO L J Ju/c 



In the Shockley-Quiesser limit [59] of ideal PV cell, the only recombination 
of the charge carriers happens through the radiative processes so the electron 
current is simply: j e = ej p h- Then the electrical power generated in the PV 
cell is Ppv = j e V , and efficiency of device is 

Ppv ejphVo 
Vtpv = p— = -5 • (6.14) 

Let us now choose the graphene's chemical potential to be \x = 0.25 eV, and 
that the PV cell is held at room temperature T2 = 300K and distance D — 10 
nm away from the graphene sheet. Then for the case of graphene's temperature 
T\ = 600 K and biased voltage V Q = 0.08 V, the output power density of our 
TPV device is Ppy/A = 6 W/cm 2 with an efficiency of rj = 35%. We note that 
these are remarkably high power densities considering that our thermal emitter 
is only one atom thick. To get a better sense of the scales involved we can 
compare the far field radiative power exchange between two black bodies 
held at temperatures T\ and T2 but involving only the photons of energies 
above the given PV band gap, and the near field radiative power exchange 
Plod between graphene and a PV cell held at these temperatures. For the 



general, the optimal voltage V depends on other parameters in the system. We avoid 
the full optimization procedure, and, motivated by the observed dependence of efficiency on 
V , choose a voltage slightly below the limit V™ ax = u g (l - T 2 /T 1 ). 
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temperatures T\ = 600 K and T 2 = 300 K one obtains P^ ad /P r ^[f = 62 times 
increase over the black body case. 

Further on, note that the near field heat transfer is particularly convenient 
since the energy is transfered by the evanescent modes and the photons with 
energy below the band gap, that are not absorbed by the PV cell, simply 
return to the graphene emitter as heat, unlike the far field case where they 
are lost in the form of propagating waves. This results in the high efficiencies 
r) = 35%, however note that these numbers are still below the Carnot limit 
7] = 50% for temperatures T\ = 600 K and T 2 = 300 K. The reason for this is 
the broad band plasmon spectrum contributing to the heat transfer with the 
high energy photons (u> > uo g ) wasting the energy difference (AE = huo — huj g ) 
on the thermalization losses heating up the system. 

To achieve even higher efficiencies one would need to tailor the emitter prop- 
erties so that it selectively radiates only in the small interval around the band 
gap of the PV cell. One way to do this would be to use surface plasmons at 
metal-dielectric surface, since they have very large density of states around the 
surface plasmon resonance; see figure l3~Tl b. and compare it to broad band spec- 
trum of graphene plasmon mode from figure l3T2l d. However, the problem with 
metals is that surface plasmon resonance usually falls in the visible/ultraviolet 
regime which is impossible to excite thermally. Alternatively one could use 
highly doped semiconductors like Indium-Tin-Oxide [70] which has a reso- 
nance in the infra-red [71], however due to high doping level there is a lot of 
electron-impurity scattering and high losses result in reduced efficiencies. In 
that regards graphene TPV system shows large promise for a new temperature 
range (600 — 1200K) solid state energy conversion, where conventional thermo- 
electrics can not operate due to high temperatures and far field TPV schemes 
suffer from low efficiency and power density. 
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Chapter 7 



Summary 



We have explored light-matter interaction in graphene in the context of plas- 
monics and other technological applications but also used graphene as a plat- 
form to explore many body physics phenomena like the interaction between 
plasmons, phonons and other elementary excitations. Plasmons and plasmon- 
phonon interaction were analyzed within self-consistent linear response ap- 
proximation. We demonstrated that electron-phonon interaction leads to large 
plasmon damping when plasmon energy exceeds that of the optical phonon but 
also a peculiar mixing of plasmon and optical phonon polarizations. Plasmon- 
phonon coupling is strongest when these two excitations have similar energy 
and momentum. We also analyzed properties of transverse electric plasmons 
in bilayer graphene. Finally we have showed that thermally excited plasmons 
strongly mediate and enhance the near field radiation transfer between two 
closely separated graphene sheets. We also demonstrated that graphene can 
be used as a thermal emitter in the near field thermophotovoltaics leading to 
large efficiencies and power densities. Near field heat transfer was analyzed 
withing the framework of fluctuational electrodynamics. 

In Chapter [2] we presented analytical methods that were used throughout 
the text. We have derived electron band structure and electron-phonon inter- 
action using the tight binding approximation. After that we derived the linear 
response functions (density-density and current-current) and used the fluc- 
tuation dissipation theorem to calculate the current-current correlation func- 
tion induced by the thermal fluctuations in the system. Finally we employed 
these results to calculate radiative heat transfer between two graphene sheets. 
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In Chapter [3] we have investigated plasmons in doped graphene and demon- 
strated that they simultaneously enable low-losses and significant wave local- 
ization for frequencies of the light smaller than the optical phonon frequency 
fouoph ~ 0.2 eV. Interband losses via emission of electron-hole pairs (1 st order 
process) were shown to be blocked by sufficiently increasing the doping level, 
which pushes the interband threshold frequency oj in te r toward higher values 
(already experimentally achieved doping levels can push it even up to near 
infrared frequencies). The plasmon decay channel via emission of an optical 
phonon together with an electron-hole pair (2 nd order process) is inactive for 
u < Uoph (due to energy conservation), however, for frequencies larger than 
Uoph this decay channel is non- negligible. This is particularly important for 
large enough doping values when the interband threshold a; inter . is above ojoph'- 
in the interval LOoph < v < Winter the 1 st order process is suppressed, but the 
phonon decay channel is open. 

In Chapter H] we showed that graphene can also support unusual transverse 
electric plasmons and we predicted the existence of TE plasmons in bilayer 
graphene. We found that their plasmonic properties are much more pro- 
nounced in bilayer than in monolayer graphene, in a sense that the wavelength 
of TE plasmons in bilayer can be smaller than in monolayer graphene at the 
same frequency. 

In Chapter we analyzed the coupling of plasmons with intrinsic optical 
phonons in graphene by using the self-consistent linear response formalism. 
We found that longitudinal plasmons (LP) couple only to transverse optical 
(TO) phonons, while transverse plasmons (TP) couple only to longitudinal 
optical (LO) phonons. The LP-TO coupling is stronger for larger concentration 
of carriers, in contrast to the TP-LO coupling (which is fairly weak). The 
former could be measured via current experimental techniques. Thus, plasmon- 
phonon resonance could serve as a magnifier for exploring the electron-phonon 
interaction in graphene. 

In Chapter [6] we analyzed the near field heat transfer between two graphene 
sheets mediated by thermally excited plasmon modes, and we demonstrated 
that there is a large enhancement of heat transfer compared to the far field 
black body radiation. Finally we showed that graphene can be used as a 
thermal emitter in the thermo-photo-voltaic system resulting in high device 
efficiencies and power densities. 
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Appendix A 

Plasmon-phonon coupling in the 
context of Feynman diagrams 



In this appendix we give alternative derivation of plasmon-phonon coupling 
in the context of Feynman diagrams. In that respect let us start by writing 
Coulomb potential 



e 2 



V(t) = - , (A.l) 

47re r 

and its Fourier transform in two dimensions 

v ™ - h (A ' 2) 

Bare Coulomb interaction V(q) can polarize electron gas by creating electron- 
hole pair which in turn screens the bare interaction resulting with an effective 
interaction W(q, u). This process can happen several times in a row (see figure 
lA.ip so we can write self-consistent equation for the effective interaction 



[-iW(q, to)] = [-iV(q)} + [-zV(q)] [-zTI(q, u>)] [-iW(q, to)] . (A.3) 

If we now use the Random Phase Approximation which neglects higher order 
scattering of the created electron-hole pair, then the polarizability Il(q, a;), 
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(b) (d) 

k + q,a>+v 



k,a> 




Figure A.l: (a) Feynman diagram for bare Coulomb interaction V(q). (b) 
Polarizability II(q, u). (c) and (d) Screened Coulomb interaction W(q,oj) in 
the Random Phase Approximation. 



depicted with a Feynman diagram in figure IATT1 (b), can be written as 

—-- ^2 iG (ri, k + q, u + u)iG (n, k, v) 

n,n' 

x (n'k + q|e iq >k)(nk|e- iqi Vk + q). (A.4) 

Here |nk), i.e. the wave function ^> n k(r) = (r|nk), denotes single particle free 
Dirac electron states (see relation (I2.29P ) and the Green function Go(n,\£,v) 
is given by expression [66] 

Go(n, k, v) = - l ~J n \. + - h , (A.5) 

fty - # nk + If] hv - E nk - IT] 

where f n y. denotes the Fermi-Dirac distribution. After integration over energy 
v we obtain 

lCl ' Wj ~ 7 (27T) 2 ^ ^ + £ nk - K'k+q 

x (n'k + q|e iqr |nk><nk|e -<qr |n'k+q), (A.6) 
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and then by using the exact wave function for Dirac electrons ?/>nk(r) = (r|nk), 
given in equation (12.291) . we obtain the polarizability 

TTY \ — A [ ^ \^ /nk ~ /n'k+q 

lq ' U) ~ (2tt)2 2^ hu + E nk - £ n , k+q 

n,n' 

x -(1 + nn' cos[^(k + q) - <p(k)]). (A.7) 
Further on, by using relation fl A. 3j) we can write the screened interaction 

where we recognize the expression for dielectric function of electron gas 

e(q, u>) = 1 - Il(q, u)V(q) = 1 - -^-1%, u). (A.9) 

Ze Q q 

Finally, we note that plasmons are simply defined as zeros of the dielectric 
function: e(q, u) = 0. 



Let us find now the phonon Green function for free phonons at zero tem- 
perature. Since longitudinal and transverse optical phonons are degenerate at 
energy huiQ = 0.196 eV, then the Green function for both branches is given by 

D °M> w ) = i Pi Pi ^2 ■ (A. 10) 

Now, the electron-phonon interaction was given in equation (12. 62p 

H e _ ph = L 2 Y,9M^pXQw (A-ll) 

The phonon motion can in turn polarize the electron gas which is described by 
a self-consistent equation for the phonon Green function renormalization (see 
figure IA~2l 



[-iD^cu)} = [-iD0(q,w)] + [-iD°(q,w)] [-zn e _ p/l (//, q, w)] [-iD M (q, w)] , 

(A.12) 



S3 



so the renormalized Green function is given by an expression: 

£>!(q,w) 



2hWo (A.13) 



(hw) 2 - (hw ) 2 - 2^a; n e _ p / 1 ( / u, q, u) ' 

Finally, the renormalized phonon frequency is defined by a pole of the Green 
function 

u 2 -u 2 = 2^-U e _ ph {fi,q,u). (A. 14) 

Up to the lowest order, the interaction will create virtual electron-hole pair 
(see figure IA.2I) which can be described with a polarization function 

— — - ^2 iG o(n', k + q, u + v)iG (n, k, v) 

n,n' 

x (n'k + qlMq^e^lnkXnklAf^e-^ln'k + q). 

(A.15) 

We note here that M qM given by equation (12.6ip is two by two matrix so that 
polarizability 11° h (fi,q, cu) isn't simply proportional to the function Il(q, ui) 
which was obtained in relation to the screened Coulomb interaction (see re- 
lation (]A.4jl ). In the context of Feynman diagrams we can say that diagram 
vertices are different for the case of Coulomb (electron-electron) interaction 
from the case of electron-phonon interaction. Further on, we obtain 



n n' 



frik /n'k+q 



tiW + — -En'k+q 

x (n'k + q|M qM e iqr |nk)(nk|M^e^ qr |n'k + q). (A. 16) 

Let us now take the exact form of the wave function ip n ^ and the matrix 
element M m according to relations (12.291) and (12.61]) . We can see that the 
polarizability depends on the phonon polarization and we obtain 



ttO it \ a 2 f d-k ^-^ /nk — /n'k+q 

p J (2vr) 2 ^ hu + E nk - £vk+ q 



1 

x -(1 - nri cos[2y?(q) - y?(k) - y?(k + q)]), and (A.17) 
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Figure A. 2: (a) Feynman diagrams for phonon Green function renormal- 
ization. (b) Feynman diagrams for polarizability function. Note that the 
electron-electron interaction vertex is different from the electron-phonon ver- 
tex. 



n e _ ph (T,c l ,u) = v 



/rak /n'k+q 



(27r) 2 ^ frw + E nk - £ n ' k+q 



x -(1 + nri cos[2<^(q) - p(k) - tp(k + q)]). (A.18) 



If we now imagine that phonon energy and momentum matches plasmon 
energy and momentum, then the electron-phonon interaction will be amplified 
through the collective electron response. In that case it won't be sufficient to 
calculate only the polarization of single electron hole pair and we will have 
to take into consideration contribution from the infinite sequence of bubble 
diagrams (which are in fact necessary to describe plasmon excitation). The 
easiest way to do this is to take the electron-phonon interaction which polarizes 
a single electron hole pair and include the possibility that Coulomb interaction 
can in turn create another electron hole pair. The infinite sequence of diagrams 
can be included if we work from the start with screened Coulomb interaction 
instead of bare interaction and one should take special account of the nature of 
diagram vertices considering if the electron hole pair was created by Coulomb 
or electron-phonon interaction (see figure IA.2I) . In that way we obtain the 
screened electron-phonon polarizability in the Random Phase Approximation 

H-e-ph(v, q, w) = n°_ ph (/i, q, u)+nl_ ph {n, q, u)W(q, u)U 2 e _ ph {fi, q, u). (A.19) 
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Here W(q,u) = x„n(q^)y(q) so we immediately see that if phonon disper- 
sion crosses the plasmon dispersion then the electron-phonon interaction will 
be amplified due the collective electron response where we have e(q, u>) = 
1 — Il(q, u)V(q) = 0. That part is in fact responsible for the plasmon phonon 
coupling. Finally the polarizability describing the bubble with different ver- 
tices is given by 

ttI ( \ tt2 / \* a [ v fnk ~ /n'k+q 

K-A^ q. ") = K e - P M q, a;) = AgJ ^ ^ ^ + — - — 

n,n' H 

x (n'k + q|M q ^e iqr |nk)(nk| e - iq V k + q). 

(A.20) 

If we include here the exact wave function i/j n k and matrix elements M qM , we 
obtain different expressions depending on the phonon polarization: 



rrl / r \ a f ^ fnk — / n 'k+q 



2 

X-(n sin[yj(q) - <p(k)] + n' sin[y?(q) - </>(k + q)]), and 

(A.21) 



TTl (T \ A f / nk /n'k+q 



1 

X-(n cos[^(q) - <f(k)} + n cos[^(q) - (p(k + q)]). 

(A.22) 

Let us first analyze interaction with the longitudinal optical phonons. In that 
respect let us take expression ( 1A.21I) and assume, without the loss of generality, 
that y?(q) = i.e. vector q is along x direction. We than obtain 



rrl / r \ A f dk v fnk /n'k+q 



(2vr) 2 ^ frw + E n]l - 4< k+q 

n,n' M 

x ^(-n sin[^(k)] - ri sin[y?(k + q)]). (A. 23) 

But the function under the integral sign is odd with respect to reflection across 
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the x axis, meaning that the entire integral vanishes i.e. H^_ ph (L, q, u) = 0. 
In other words we have shown analytically that there is no whatsoever cou- 
pling of plasmons and longitudinal optical phonons! Finally, to find the cou- 
pling of plasmons with transverse optical phonons one only needs to solve self- 
consistent set of equations ()A.14j) and (IA.19[) which was done numerically and 
demonstrated to agree with the results of chapter [5j where we used different 
gauge to obtain the same result. 
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List of Figures 



2.1 a) Graphene crystal structure. Unit cell vectors are a x and a 2 
while A and B are atoms of the basis, b) Brillouin zone. We 
mark high symmetry points K and K' where the low-energy 
electron excitations are described by massless Dirac equation. 

2.2 Graphene electron band structure with Dirac cones around K 
point (magnified). Intrinsic graphene has Fermi level Ep = 
E K = 

2.3 a) Atom motion in graphene during the longitudinal optical 
phonon oscillations, b) Motion of basis atoms described by a 
vector u in the real space induces vector potential A that moves 
Dirac points in the reciprocal space. Dirac point symmetry 
causes unusual polarization of this vector potential: Alu. 

2.4 (a) Schematic diagram of the radiation transfer problem: a free 
standing sheet of graphene at temperature T\ is radiating to 
another free standing graphene sheet at temperature T 2 and 
distance D away, (b) Polarization vectors defined in the text. 

3.1 (a) Schematic description of a surface plasmon (SP) on metal- 
dielectric interface, (b) SP dispersion curve (solid blue line) 
for Ag-Si interfaces; dotted blue is the light line in Si; dashed 
red line denotes the SP resonance, (c) Wave localization and 
propagation length for SPs at Ag-Si interface (experimental Ag 
losses are taken into account) 
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2 (a) Schematic of the graphene system and TM plasmon modes. 
Note that the profile of the fields looks the same as the fields of 
an SP [Fig. 13.1( a)], (b) Electronic band structure of graphene; 
to indicate the vertical scale we show the Fermi energy level 
for the case Ep = 1 eV. (c) Sketch of the intraband (green 
arrows) and interband (red arrows) single particle excitations 
that can lead to large losses; these losses can be avoided by im- 
plementing a sufficiently high doping, (d) Plasmon RPA and 
semiclassical dispersion curves. Black solid (RPA) and black 
dot-dashed (semiclassical) lines correspond to e r i = e r 2 = 1; 
Blue dashed (RPA) and blue dotted (semiclassical) lines cor- 
respond to e r i = 4 and e r 2 = 1. The green (lower) and rose 
(upper) shaded areas represent regimes of intraband and inter- 
band excitations, respectively 

3 Properties of plasmons in doped graphene. Solid-lines are ob- 
tained with the number-conserving RPA calculation, and the 
dashed lines with the semiclassical approach. Losses (a), field 
localization (wave "shrinkage") (b), and group velocity (c) for 
doping Ep = 0.135 eV, and relaxation time r = 1.35 x 10~ 13 s, 
which corresponds to the mobility of 10000 cm 2 /Vs. The up- 
per scale in all figures is frequency v = u/2tt, whereas the rose 
shaded denote the region of high interband losses. . . . 

4 (a) The real part of the conductivity in units of <Jq = ne 2 /2h 
in dependence of frequency fiw/Ep, and (b) the corresponding 
relaxation time as a function of wavelength. The contribution to 
$t,a(uj) from impurities is chosen to be negligible. The displayed 
graphs correspond to two different values of doping which yield 
Ep = 0.135 eV (solid blue line), and E F = 0.640 eV (dashed red 
line). The position of the optical phonon frequency fajJoph ~ 0.2 
eV is depicted by the dotted vertical line in (b); dot-dashed lines 
depict the values of wavelengths corresponding to 2Ep, that is, 
the interband threshold value (for q = 0) for the two doping 
concentrations 
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3.5 Properties of plasmons in doped graphene. Solid-lines are ob- 
tained with the number-conserving RPA calculation, and the 
dashed lines with the semiclassical approach. Losses (a), field 
localization (wave "shrinkage") (b), and group velocity (c) for 
doping E F = 0.64 eV; losses are calculated by using the relax- 
ation time = r D ^ + r~_ h , where tdc = 6.4 x 1CT 13 s, and 
T e -ph is the relaxation time from the electron-phonon coupling 
for the given parameters. In the white regions (right regions in 
all panels), losses are determined by tdc- In the yellow shaded 
regions (central regions in all panels), losses are determined by 
the optical phonon emission, i.e., r e - p h- The rose shaded areas 
(left region in all panels) denote the region of high interband 
losses. Dotted vertical lines correspond to the optical phonon 
frequency uo P h ~ 0.2 eV. The upper scale in all figures is fre- 
quency v = u/2tt. See text for details |49 



4.1 The band-structure of bilayer graphene. The two upper bands 
(as well as the two lower bands) are perfectly nested and sepa- 
rated by 7 ~ 0.4 eV; q = r y/hv F . Horizontal line depicts one 
possible value of the Fermi level, and arrows denote some of the 
possible interband electronic transitions. See text for details. . [54 



4.2 The real (red dotted lines) and imaginary (blue solid lines) part 
of the conductivity of bilayer graphene for two values of doping: 
fx = O.47 (a), and \i = O.97 (b). The conductivity is in units 
of cr = e 2 /2h, and the frequency is in units of u = ■y/k. The 
5-functions in 3ftcr(u;) at oj = (intraband transitions) and u = 
~f/h (transitions from the lower to the upper conduction band 
depicted as green solid arrows in Fig. 14.11) are not shown (see 

!) E 
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4.3 The plasmon dispersion curve Aq = q — oj/c vs. u for fi = 
O.47 (a), and /i = O.97 (b) is shown as blue solid line. To the 
right of the vertical red dotted lines plasmons can be damped 
via excitation of electron-hole pairs, whereas to the left of this 
line these excitations are forbidden due to the Pauli principle. 
Black dashed line in (b) (which closely follows the blue line) 
corresponds to Eq. (j4.8j) . The wave vector is in units of go = 
■y/hvp, and the frequency is in units of u = j/h 



5.1 (a) Schematic illustration of the lattice structure with two sub- 
lattices (A and B). (b) The two degenerate Dirac cones are 
centered at K and K' points at the edge of the Brillouin zone. 

(c) A displacement of lattice atoms u(r) is parallel (perpendic- 
ular) to the propagation wave vector q of a LO (TO) phonon. 

(d) The displacement u(r) creates an effective vector potential 
A c g perpendicular to u(r) (the sign of A e g for the K' point is 
opposite to that for the K point) 



5.2 Dispersion lines of hybrid LP-TO plasmon-phonon modes (solid 
lines) and of the uncoupled modes (dashed lines) for two values 
of doping: (a) n = 5 x 10 12 cm -2 , and (b) n = 5 x 10 13 cm -2 . 
The hybridization is stronger for larger doping values. Grey 
areas denote the region of single-particle damping 



5.3 Dispersion lines of hybrid TP-LO plasmon-phonon modes (solid 
lines) and of the uncoupled modes (dashed lines) for two values 
of doping: (a) n = 7.5 x 10 11 cm -2 , and (b) n = 9.5 x 10 11 cm -2 . 
The plasmon-like dispersion is very close to the light line q = 
oj/c; therefore, the ordinate shows Aq = q — oj/c 



98 



6.1 (a) Schematic diagram of the radiation transfer problem: a free 
standing sheet of graphene at temperature Ti is radiating to an- 
other free standing graphene sheet at temperature T 2 and dis- 
tance D away, (b) Schematic diagram of the field profile for even 
mode, (c) Odd mode, (d) Contour plot of the transfer function 
h^f for the case of two graphene sheets at the same chemical 
potentials /ii i2 = 0.5eV and same temperatures Ti |2 = 300A', 
separated by a distance of D = 10 nm. Dashed line denotes the 
plasmon dispersion relation for a single isolated graphene sheet, 
while poles of transfer function h p n ^ show dispersion relations of 

the coupled (even and odd) modes of the two sheets |72 

6.2 Contour plot of the near field heat transfer between two graphene 
sheets H™f normalized to the far field heat transfer between two 
black bodies H% B of the same temperatures, in the log scale. 
Here T 2 = 300A, // 1>2 = O.leV and r 1>2 = 10~ 13 s [75 

A.l (a) Feynman diagram for bare Coulomb interaction V(q). (b) 
Polarizability II(q, u). (c) and (d) Screened Coulomb interac- 
tion W(q, to) in the Random Phase Approximation |82 

A. 2 (a) Feynman diagrams for phonon Green function renormaliza- 
tion. (b) Feynman diagrams for polarizability function. Note 
that the electron-electron interaction vertex is different from the 
electron-phonon vertex |85 
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